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SECTION 1 


INTRODUCTION 


Volume 2 of the present report giver an account cf theoretical 
advances and associated demonstrations made during the reporting period 
relating to the synthesis of active controllers for large flexible 
structures. Particular applications discussed include the quenching of 
vibrations in support structures for optical pointing systems and the 
efficient generation of closed-loop large-angle slewing maneuvers for 
linear tracking. Certain aspects of discrete-time dynamic modeling for 
structural vibration control are also treated in depth. The vehicles for 
the demonstration examples are the ACOSS Model No. 2 representation of an 
optical pointing system (for vibration control) and a rotating rigid hub 
with four flexible appendages (for large-angle slew control). A concise 
summary of the principal results contained in the present volume is given 
in the remainder of Section 1. 

1.1.1 Vibration Control 

Effective accommodation of disturbances is fundamental to closed- 
Ictp regulation, being especially vital for the quenching of vibrations 
in large flexible optical pointing systems. Accommodation of a (small) 
finita number of periodic disturbances at fixed known frequencies poses 
only mild difficulty. Substantially more difficult is the problem of 
accommodating an aperiodic disturbance having nonnegligible power 
spectral density across a broad range of frequencies. In the preceding 
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reporting period, a scientific experiment was conceived to examine the 
question of active control synthesis for broadband disturbance accommoda¬ 
tion in which the selection of reduced-order models and of active 
transducers were to be given substantial consideration. 

Sections 2 and 3 of the present volume taken together constitute 
an initial account of progress in the study of broadband disturbance 
accommodation for large flexible optical pointing systems. The 
disturbance to be neutralized has a constant nonzero power spectral 
density over the frequency range between 0 and 5 Hz. In the first part 
(Section 2), attention is focused primarily upon the selection of active 
transducers, both by functional type and by location. The theoretical 
basis of a recently-developed selection algorithm is fully described; the 
principal elements are full-rank matrix factorizations -- in particular, 
the QR-decomposition — and their employment in the representation of 
solutions to least-squares approximation problems. The algorithm is 
designed to select a relatively small number of actuators that have a 
relatively high influence on the variables to be regulated (in this case, 
optical system pointing errors). Results of four examples with ACOSS 
■todel No. 2 demonstrate that selections are readily achievable that are 
substantially fewer in number than the number of controlled modes, and 
yet give promise of effective performance. In the second part (Section 
3), attention is shifted to the design and evaluation of controllers for 
ACOSS Model No. 2 associated with the four active transducer selections 
described in Section 2. The controller feedback structure is generated 
using a textbook disturbance-rejection control design. Attempts at 
stability-enhancing adjustments beyond meeting the specifications on the 
optical pointing error are deliberately postponed in order to isolate the 
influence upon the overall synthesis-process of the selection of reduced- 
order models and of active transducers, including their mutual inter¬ 
action. In each example, the effect of unmodeled modes upon control 
system scability and performance is examined using an expanding family of 
reduced-order evaluation models. In the evaluations, closed-loop 
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stability over a range of frequencies substantially beyond that of the 
design model is demonstrated. 

In the modeling of dynamic systems, state space representation 
provides a unified framework which abstracts certain system-theoretic 
characteristics common to a wide variety of physical systems. The 
resulting separation of dynamic characteristics from specific physical 
processes provides beneficial clarification. However, it has the 
disadvantage of masking the physics of the problem. The sharpened 
understanding that follows from taking into account the unique 
characteristics of a specific physical process usually increases the 
performance achievable from a process controller. In this spirit, a 
frrsh approach to sampled-data control for flexible structures* is taken 
in Section 4 i~ which second-order difference equations provide the 
dynamic model. For the problem of velocity output feedback, explicit 
sufficient conditions for stability in terms of the feedback gain matrix 
and the sampling rate are given. Surprisingly, these conditions require 
non-colocation of actuators and sensors as a function of the sampling 
rate. The non-colocation property collapses to the familiar colocation 
property for analog systems as the sampling period converges to zero. 

1.1.2 Large-Angle Slew Control 

Large-angle slewing maneuvers for flexible spacecraft generally 
require the consideration of nonlinear structural dynamics. However 
certain problems of substantial interest for applications require only 
the consideration of linear autonomous structural dynamics. Several such 
pfotr.ems are formulated and solved in Section 5; namely, the problem of 
optimally tracking a known trajectory (without hard terminal 
constraints), and the problem of optimal linear slew in which the 
terminal values of some or all of the state variables at a fixed finite 
terminal time are prescribed. The emphasis is upon obtaining explicit 
expressions for solutions of the necessary conditions for optimality in 
terms of matrix parameters which are solutions of linear algebraic 
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equations. Hie problem of computing solutions is thereby reduced from 
one of numerical calculus (i.e., integration) to one of numerical linear 
algebra. Significant increases in computational efficiency resulting 
from this approach are demonstrated with a variety of examples on a 
rotating rigid hub with four flexible appendages. The method of 
obtaining explicit solution expressions is also shown to have direct 
application to the efficient generation of feedback control solutions to 
the general unconstrained fixed time linear optimal control problem. 

1.2 Limitations 

The scientific experiment relating to active control synthesis for 
broadband disturbance ac<- .nmodation is incomplete. Projected work 
includes: 

(1) Expansion of the class of candidates from which active 
transducers are selected; 

(2) Examination of the interaction between the selection of 
reduced-order models and of active transducers; and 

(3) Modifications of the controller feedback structure to ensure 
an acceptable compromise between the requirement to 
neutralize the effects of disturbances, and the limitations 
imposed by an incomplete knowledge of structural and 
disturbance characteristics. 









PROGRESS ON SYNTHESIS OF ACTIVE CONTROL FOR BROADBAND 


DISTURBANCE ACCOMMODATION 

PART 1: ACTIVE TRANSDUCER SELECTION 


2.1 Background 

In the automatic control of dynamic systems, the achievement of 
stability and a specified Level of performance, in spite of anticipated 
disturbances, is fundamental. When the object of control has substantial 
inherent flexibility, the problem of accommodating disturbances is 
aggravated by the tendency of the disturbances to excite structural 
vibrations. Even with flexible structures, the effects of.a few periodic 
disturbances at known frequencies are relatively easy to accommodate. 
However, when the frequency of periodic disturbances is unknown or 
variable, or when the disturbances are aperiodic, disturbance 
accommodation i3 much more difficult. Relatively little has been written 
regarding accommodation of the latter class of disturbances in the 
context of large flexible structures. One careful study treating 
unknown-frequency periodic disturbances [Ref. 2—1j has reported moderate 
success, l't is notable that the treatment reported gave careful 
attention to the placement of active transducers (i.e., actuators and 
sensors). The only study available to date on accommodation of aperiodic 
disturbances with a continuum frequency spectrum ("broadband") in the 
context of large flexible structures [Ref. 2-2) also reports moderate 
success, without however, giving any attention to the selection of active 




transducers beyond the analysis developed for accommodating periodic 
disturbances of known frequencies [Ref. 2-3]. 

The presence of a broadband disturbance escalates the importance 
of giving careful attention to the selection of reduced-order models, the 
selection of active transducers, and the interdependence of tnese two. 
selections. In Reference 2-4, we outlined a scientific experiment (not 
in hardware) for the purpose of systematically investigating the impact 
upon the performance of active broadband disturbance accommodation of 
various assunptions relating to the selection of reduced-order models and 
of active transducers. This ex eriment is not a single design attempt, 
but rather a process involving a sequence of many variously-constrained 
designs, in order to converge on an appropriate overall algorithm for 
broadband disturbance accommodation with large flexible structures. 

Sections 2 and 3 of the present report constitute the initial 
report of progress in the conduct of the aforementioned experiment. 
Section 2 focuses principally on the theory and application of the active 
transducer selection strategy being employed, while Section 3 focuses on 
the related matters of controller design and evaluation. Two very 
encouraging trends have emerged from the results reported herein: 

(1) Under appropriate assumptions, substantial reductions in the 
number of active transducer pairs below the number of controlled modes 
are feasible, which implies substantial design flexibility may be 
available; and 

(2) Using an expanding family of evaluation models, stability of 
reduced-order controllers over a frequency range well beyond the range of 
the design model has been demonstrated. 

The significance of item (2) is that sophisticated adjustments of the 
controller design have been deliberately postponed; this result reflects 
principally the reward of careful attention to active transducer 
selection. 
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The remainder of Section 2 is briefly summarized. The active 
transducer selection strategy being employed in the scientific experiment 
was briefly outlined in a previous report [Ref. 2-5, Sec. 3]. In Section 
2.2, a complete account of the theory that underlies this strategy is 
given. In addition, all phases of the strategy are fully discussed in 
the light of the theoretical building blocks. In Section 2,3, a number 
of applications of the selection strategy to ACOSS Model No. 2,.under 
varying assumptions, is presented. In Section 2.4, some principal 
observations on the results are made. An appendix, Section 2.5, contains 
the proof of some technical matters that arise in Section 2.2. 


2.2 Theory: A Selection Algorithm 

Most strategies for active transducer selection reported in the 
literature use appropriately-defined indices of controllability and 
observability to evaluate alternatives. The selection strategy proposed 
and briefly outlined previously [Ref. 2-5, Sec. 3] approaches the 
selection problem from a different direction. Linear combinations of the 
modal influence vectors associated with each actuator are sought which 
best approximate the modal coefficient vectors of variables to be 
regulated (e.g., line-of-sight errors in an optical system). The 
principal theoretical building blocks of this strategy are the special 
properties of full-rank matrix factorizations (in particular, the 
QR-decomposition) and the representation theory for solutions of least 
squares approximation problems. In what follows, a complete though 
concise account of this underlying theory is necessary to enable a proper 
description of the construction and operation of the selection algorithm. 






2.2.1 Fundamental Concepts 

As will be seen shortly, least-squares approximation problems lead 
naturally and directly to the Moore-Penrose generalized inverse. This 
construction exists for arbitrary rectangular matrices A: nxm, but is not 
quite so direct in the case that the inequality 


k = rank(A) £ min {n,m} (2-1) 

becomes strict,. In that case, the ability to decompose A into factors is 
extremely helpful. A decomposition of the form 

A = BC (2-2) 


where E: nxk and C: kxm satisfy 

rank(B) ■ k » rank(C) 

with k defined by Eq. (2-1), is called a full-rank factorization of A. 
It is well known that every matrix has such a factorization [Ref. 2-6]. 
Although such factorizations are not unique, the class of all such 
factorizations is generated from a known factorization by the relations 

B 1 =» BY , » Y _1 C 


for some nonsingular matrix Y [Ref. 2-7]. Full-rank factorizations 
assune great importance in numerical linear algebra [Ref. 2-8]. 

Essential to the selection algorithm to be described later is the 0R- 
decomposition [Ref. 2-9], which always has embedded within it a full-rank 
factorization. 
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2.2.1.1 Ohe QR-Decomposition 

The classical Gram-Schmidt orthcgonalization [Ref. 2-10, pp. 

127—1261 is a constructive process of transforming a set of linearly 
independent vectors in an inner-product space into an orthonormal set. 
The matrix representation of the result is the QR-decomposition. It 
follows that every matrix has a QR-decomposition, whose specific 
structure depends on the rank of the matrix. Since the results are 
thoroughly discussed in the literature [e.g.. Ref. 2-11], the parts 
relevant here are simply quoted. 

Theorem 2-1 (Existence of QR-Decomposition). Assume that the matrix A: 
nxm has rank k=min{n,m}. Then there exists a (nonsingular) column 
permutation matrix P: mxm, a unitary matrix Q: nxn, and a matrix R: nxm 
containing an upper triangular submatrix R^: kxk with positive diagonal 
elements, such that 


The structure of the R matrix is 


[*.] 



(2—4a) 


or 

[*, : s ] 


(2-4c) 


according as A is determinate (k=n=m), overdetermined (k=m<n), or under¬ 
determined (k=n<m), respectively. □ 
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The word "unitary" refers, in the complex case, to a matrix allowing 
complex entries whose columns are orthogonal relative to the complex 
inner product; the corresponding property in the reaj- case refers to a 
matrix restricted to real entries whose columns are orthogonal relative 
to the real inner product ("orthogonal"). 

, Theorem 2-1 has an obvious extension to the rank-deficient case 
that k < minfn,m). The extended result is used only in the proof of 
Proposition 2-8 in the Appendix (Section 2.5). The decomposition of Eq. 
(2-3) is a full-rank factorization (cf. Eq. (2-2)) except in the over¬ 
determined case (cf. Eq. (2-4b)). In that case, Eq. (2-3) is equivalent 
to the full-rank factorization 

AP =» QiRi 

where Qi is the partition of Q compatible with R-| . 

For the intended application, stronger properties of the 
nonsingular submatrix R-j are required than are evident in the 
Gram-Schmidt interpretation as reflected in Theorem 2-1. A powerful tool 
for accomplishing this is the following triangularization lemma of 
Householder [Ref. 2-12]. 


Lemma 2-2 . Assume that n-vectors a*0 and v are specified. Then there 
exists a unit vector w: nxl such that 

H(a) = (I - 2 ww*)a = lalv □ (2-5) 

n 

In Eq. (2-5), I n is the nxn identity matrix, the asterisk (*) denotes 
the complex conjugate transpose operation, and the double bars denote the 
vector norm generated by the inner product: 
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(a,a) 


a a 


Successive application of Householder transformations, each having the 
form defined in Eq. (2-5), to columns of the A matrix, including 
appropriate rearrangement (pivoting) of the columns at each step, leads 
to the desired triangularization. Reliable and widely available software 
for computing a QR-decomposition uses the Householder lemma [Ref. 2-13J. 
The stronger property that results is the following. 


Theorem 2-3 . In Theorem 2-1, elements of the matrix Ry = [rj,j ] may 
be taken to have the diagonal dominance property : 


> Id, 




(2-6) 


In particular, this implies the (weaker) property that the diagonal 
elements form a nonincreasing sequence: . 


U - 1+ 1,A+1 


The subtle distinction between the two theorems above should 
be carefully noted. The results of Theorem 2-1 dascribe every 
QR-decomposition; the sharper results of Theorem 2-3 describe only those 
QR-decompositions generated by the use of Householder transformations. 

2.2.1.2 Solutions to Least-Squares Problems 

Least-squares problems for finite-dimensional linear systems of 
the form 
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(2-8) 


Ax - b 

with A: nxm become significant in the overdetermined case (n>m) in which 
there are too many equations for a solution to exist (in general). 
Instead, one attempts to minimize a positive-definite quadratic function 
of the difference, such as 


J(x> A, b) = (b-Ax)*(b-Ax) (2-9) 

or, more generally (apparently), 

J (x; A, b) = (b-Ax)*W(b-Ax) (2-10) 

W , 

where Ws nxn is Hermitian (i.e., W*=W) and positive definite. Being 
positive definite however, W can be written in the factored form 

W ■ L*L 

with L: nxn nonsingular (Ref. 2-10, pp. 139-140], This enables Eq. 

(2-10) to be written as 

J w (x; A, b) =* (b - Ax)*(b - Ax) = J(x? A, b) (2-11) 

where b » Lb, A = LA, so that the least-squares problems associated with 
Eq. (2-10) and Eq. (2-9), respectively, are in fact equivalent under a 
simple scaling. It is therefore sufficient to discuss the relevant parts 
of the least-squares theory using Eq,» (2-9), even though the selection 
algorithm is formulated \sing Eq. (2-10). Only the unconstrained 
minimization of Eq. (2-9) is discussed. Full details, including 
treatment of constrained minimization, can be found in Reference 2-11. 

The minima of Eq, (2-9) are precisely the set of solutions to the 
normal equations 
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A (b-Ax) 


(2-12) 


It is worth noting that this result is a consequence of the geometric 
structure of inner-product spaces [Ref. 2-7] and does not follow from the 
methods of calculus (the operation x x of. complex conjugation being 
nowhere differentiable). In addition, Eq. (2-12) provides a useful 
geometric insight; at any least-squares minimum, the residual vector 

P = b - Ax 

associated with the linear system Eq. (2-8) is orthogonal to the range 
space (i.e., the subspace spanned by the columns) of A, denoted by 
col(A). 

The solutions of the normal equations form an affine space of 
dimension m-k, where m and k are defined by Eq. (2-1). Written 
explicitly, they are; 

x = A f b + (1^ - A + A)o> , w; nxl arbitrary (2-13) 

where the superscript (t) represents the Moore-Penrose inverse [Refs. 
2-14, 2-15]. The basic properties oS this generalized inverse will be 
used without elaboration, as they are thoroughly documented [Ref. 2-16] 
and have been discussed previously in a related context [Refs. 

2-17, 2-18]. It follows from these properties that the two terms in Eq. 
(2-13) are orthogonal and, therefore, that the particular least-squares 
solution 


has (uniquely) the smallest norm of all the solutions Eq. (2-13). 
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Each least-squares solution gives the minimum value 


J min “ I(I n ” AA f )bH 2 (2-14) 

for the functional Eq. (2-9). Geometrically, this represents the length 
(squared) of the projection of b onto col(A) J \ the orthogonal complement 
of the range space of A. In particular, it follows that the minimum 
valuQ is zero if and only if 


b e col(A) (2-15) 

i.e., that Eq 0 (2-8) has an exact solution. 

r ..he representation of least-squares solutions can be sharpened 
substantially if the Moore-Penrose inverse is displayed using the QR- 
Jecomposition of A. Two cases are distinguished: 

(I) k » n £ m, for which Eq. (2-15) is true; and 
(II) k = m < n, for which Eq. (2-15)'is false. 

Theorem 2-4 (Case I) . Assume that A has the decomposition of Eq. (2-3) 
with k = n <_ m. Then the minima of Eq. (2-9) are given by: 

(P _1 x) = R + Q b + (1^ - R + R)u , u: nxl arbitrary (2-16) 

If n=m, the set of Eq. (2-16) collapses to the unique solution 

(p-’x) = R’Yb • (2-17) 

In both instances, the minimum value in Eq. (2-14) is zero. □ 

This result is used in the initial stages of the selection algorithm. 
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mi even richer structure develops in the next case. 


Theorem 2-5 (Case II) . Assume that A has the decomposition of Eq. (2-3) 
with k = m < n. Denote by Q = [Q-j : Q 2 ] the partition of Q compatible 
with the partitior of R in Eq. (2-4b). Then the minimum of Eq. (2-9) is 
given (uniquely) by: 


, -1 . -1 * 

(P x) - fi, b 

and the minimum value in Eq. (2-14) is 

J *i„ ' ’Bp* 6 ' 2 ° 


(2-18) 


The significance of Theorem 2-5 arises from the riqh structure of the 
mappings associated with the matrices Q-js nxk and Q 2 : nx(n-k), which 
are briefly summarized next. They all follow from viewing Eq. (2-3) in 
the equivalent form 


Q*(AP) = R (2-20) 

which exhibits the final result of the Householder construction (cf. 
Lemma 2-2). It follows immediately that Eq. (2-20) embodies the two 
equations: 


Q/(AP) - R 1 

and 


( 2 - 21 ) 
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Q 2 *(AP) = 

0 

(2-22) 

From Eqs. (2-21) and (2-22) it follows quite directly that Q-| 

are orthogonal: 

and Q 2 


* 


«1 Q 2 " ° ' 

Q 2 Qi = 0 

(2-23a) 




«! Q 1 - ^ ' 

% °2 " X n-k 

(2-23b) 


and have range spaces 

col(Q 1 ) - col(AP), col(Q 2 ) => col( AP)"*" (2-24) 

respectively. Furthermore the product mappings 

Q Q * = I n - Q 2 Q 2 * (2-25) 

and 

Q 2 Q 2 * ■ I„ ~ Q 1 Q.,* (2-26) 


are orthogonal projections, respectively (of. Bqs. (2-23), onto the range 
spaces listed in Eq. (2-24). Note that Eq. (2-26; provides the link 
between Eq. (2-14) and the result in Eq. (2-19). Theorem 2-5 forms the 
principal basis for the selection algorithm. 

The behavior of the minimum value in Eq. (2-14) as a function of 
A, and in particular as columns are removed from A, is of considerable 
interest for the active transducer selection problem. Roughly, one 
expects that the minimum value does not change it columns of A are simply 
rearranged, and that it does not decrease if a column of A is removed. 
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Theorem 2-6 . Let p be an integer, 1 <_ p _< m, and denote by 
A il •»ip an arbitrary p-column submatrix of A. Then 

101,1 J'x; A, h) £ " an J(y; }' } , b) (2-27) 

xeE m ytsP S 

w v 3 re E a denotes the appropriate a-dimensional i n .vstt product ppace.H 


The proof of this result is deferred to the Appendix (Section 2.5). 

2,2,2 Algorithm Formulation 

It is now possible to give a precise formulation of the active 
transducer selection algorithm. It is assumed in what follows that: 

(1) A reduced-order structural model has been chosen; and 

(2) Sensors are to be colocated with actuators. 

The focus is on determining an appropriate, and relatively small, 
selection of actuators from a prespecified collection of candidates. 

Dynamic equations in modal coordinates for the reduced-order model 
to be considered have the general form: 

n + 2znn + n 2 n = (** B A ) U A + (* T V U D (2-28) 


y 


/.T_ TJ . T T, t • 

(* Cp ) n + (* ~ v ) n 


(2-29) 


(♦V' T - 


) n 


(2-30) 
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where n = (ni,•••/hn) T * s the vector of modal coordinates retained 
in the reduced-order model? fl = diag(coi): nxn is the matrix of 
characteristic frequencies; $: vxn is the truncation of the 
principal-axis matrix of transformation defining the modal coordinates 
(mode shapes of the reduced-order model), v being the number of physical 
generalized coordinates; Z = diag(^, • • »,£ n ): nxn is the matrix of 
assuned modal damping ratios representing inherent structural damping; 
u A = (ui,«*» / u m ) T is the vector of inputs to the actuators; 
up = (wi,«»»,Wy) t is the vector of disturbance inputs; 
y 2 (y, f ...,y A )T j, s the vector of outputs from the sensors; 
z = (zi,»«»,z g ) T is the vector of variables to be regulated; and 
B A ; vxm, Bq: vxy, Cp: fxv, Cy: fxv, D: 3x\> are the influence 
matrices associated with the actuators, disturbances, displacement 
sensors, rate sensors, and regulated variables, respectively. The 
superscript (T) denotes matrix transpose. Equations (2-28) through 
(2-30) are written in a form that highlights the presence of generalized 
node shapes (modal influence vectors) of the form $ T x. 

The selection problem is defined in terms of two sets of 
parameters: 

(1) - A collection of actuator candidates, consisting of the 

T i i 

generalized node shapes $ b^, where b" is the j-th column of 
the matrix B A in Eq. (2-28); and 

(2) A collection SB of regulation target vectors, consisting of 
tne generalized node shapes $ T d^, where d/ 1 is the i-th row 
of the matrix D in Eq. (2-30), 

It is assumed that the designer has a relatively wide range of choices 
for actuators (in particular, having redundancy) both by functional type 
and location for inclusion in the candidate collection^/. The idea 
behind collection SB is as follows: actuators are to be selected from 
the collection «// which have a relatively strong influence on each of the 
vectors in collection SB , thereby enhancing the prospects of achieving 
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the desired regulation of the vector z in Eq. (2-30). In order to 
clarify ideas and simplify the notation, the remaining discussion 
identifies the regulated variable z with the two-vector 
(Elqsx» e L0SY) that is of interest in the applications' of Section 2.3 
and which represents line-of-sight (LOS) error parameters for an optical 
system. The collection .35 then consists of the two elements: 


, LOSX A 
b = 


*\osx 


LOSY A 
b = 


(2-31) 


where d£ osx and d^ QSY are the rows of D corresponding to the regulated 
variables e L osx and e L0SY in (2-30). 

Using the parameters from collections ^ and .55 , an ideal 
selection problem, and a least-squares approximation to it, are defined. 
Denote by 


A = [ysJ] = [a 1 • ••• | a w ] : n x p (2-32) 

the matrix formed from elements of the actuator candidate collection . 
(The type of correspondence defined by Eq. (2-32) is used repeatedly in 
what follows.) It is assumed that the collection is rich enough to 
span the space E n of modal coordinates; i.e., 

sp(^) 5 col(A) = E n (2-33) 

Then there exist coefficient vectors x = (x .•••,x ) T , y = (y»#»,y ) T 

In In 

such that 


and 


= y U 

£j=1 


*1 a S 
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LOSY 


Ideal Selection Problem ; Find a minimal subcollection «V of and 
real coefficient vectors x' = (Xj), y' = (y j ), j e HjJ'), such that 


and 


a x 'j 


A 'y’ = Ij 


a5y ’j 


where the index set K j&') is the subset of {1 r •••#U> that identifies 
which elements of jt ara retained in vtf', and A' = [,*/']. 

As a mor? feasible alternative to the ideal selection problem, the 
following least-squares approximation is formulated.. 


Approximate Selection Problem ; Find a subcollection id' of and 
corresponding coefficient vectors x' = (Xj), y' 5 (yj), j e I(v«/') 
such that 


(1) J w (x'; A', b L0SX ) £ (b L0SX - A'x') T W(b L0SX - A'x') and (2-34) 
J H (y, A', b L0SY ) £ (b L0SY - A’y' ) T W(b L0SY - A-y) (2-35) 

are minimized, where A' = [jl' ]; and 


(2) the subcollection c zt' is minimal with respect to acceptable 
increases in the minimum values: 





j b LOSX A min a', b L0SX ) 

min x' W 


(2-36) 


(2-37) 


In the above, the matrix W: nxn is chosen to be real, symmetric, and 
positive definite. However, in view of the equivalence, under 
appropriate scaling, conveyed by Eq. (2-11), the discussion below omits 
explicit reference to the weighting matrix W. Note that, although the 
correspondence 


[j/'] = A' 

defined by Eq. (2-32) does not distinguish between distinct column 
rearrangements, Proposition 2-7 (Section 2.5) shows that such 
distinguishability is inessential for its intended usage. Thus, the 
approximate selection problem is well-defined. 

A brief overview of the algorithm for the solution of the 
approximate selection problem is as follows. A sequence of eliminations 
from the actuator candidate collection jl is generated on the basis of 
information from least-squares solutions corresponding to Eqs. (2-34) and 
(2-35) at each step. These solutions are displayed using the appropriate 
QR-decomposition. Note that each elimination also leads to a reduction 
in dimension of the inner-product space in which the least-squares 
minimization is done. No matter what elimination sequence develops, 
Theorem 2-6 guarantees a sequence of least-squares approximations that is 
steadily deteriorating in quality in the sense that the sequence of norms 
for the residual vector at each least-squares solution is nondecreasing. 
When has oeer. reduced by elimination to a collection for which 

further elimination (say to cs/") leads to the condition that either of 
the minimum values Eq. (2-36) or (2-37) (corresponding to «.*/" rather than 
to ca/' ) exceeds a specified tolerance tq? the elimination sequence is 
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terminated, and is taken to be a minimal subcollection. Next, 
details of the elimination sequence are outlined. 

2.2.2.1 Criteria for Elimination 

In view of Theorem 2-4, the search for a minimal subco2 lection in 
the sense of the approximate selection problem eventually requires 
criteria for intelligently eliminating a column from a determinate or 
overdetermined matrix. Two such criteria are discussed prior to 
describing the elimination sequence. Consider a stage of the sequence at 
which the matrix A = of the remaining actuator candidates is 

(say) overdeternuned with l < n columns. The QR-decomposition of A has 
the form of Eq. (2-3) with R having the form of Eq. (2-4b). Denote by x 
and y, respectively the unique solutions (cf. Eq. (2-18)) to the 
least-squares problems associated with Eqs. (2-34) and (2-35). The 
criteria for elimination of actuator candidates from A are defined in 
terms of the matrices R, x, and y; in particulars 

(I) the relative magnitude of the smallest diagonal element r^ 
of R^j and 

(II) the relative magnitude of the products t^Xi, i-1,•••,!, 
andrjjyj, j=1,• ••,£, respectively. 

The rationale for each criterion is briefly discussed. 


Criterion I . The scalar equations for the QR-decomposition of A have the 
triangular structure: 

* j - (2 - 38) 

where the diagonal elements rjj of Ri form a nonincreasing sequence 

(cf. Eq. (2-7)). Denote by the (£-1)-dimensional subspace spanned 

by the vectors {q 1 ,***, q*" 1 }. Since the vector a*-a*r is m S ,, 

■ ££ £-1 

is an upper bound for the Euclidean distance between a^ and S^s 
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d ( a *' s jl _ 1 ) = inf {la^-s J: s e } 


(2-39) 


Because of the triangular structure of Eas. (2-38), Sg_i is also the 
linear span of {a 1 ,***, a* -1 }. Therefore, Eq. (2-39) is a measure of 
the near-dependence of a& upon the preceding column vectors ai, 
i=1 »• •»,JM • Sufficiently small values of relative to 
suggest that a^ might be eliminated without an unacceptable increase in 
the resulting minimum value Eq. (2-36) and (2-37) associated with the 
truncated matrix. 

Criterion 11 . Eliminations in accordance with Criterion I are made only 
on the basis of the QR-decomposition of A, and without regard for the 
least-squares solutions that determine the quality of the desired 
approximation. Criterion II takes those solutions into account. In 
comparing the numbers |rj.iXj_|, i=1,•••,£, suppose that the minimum 
occurs at some index o < .1; i.e. , 


min 
i=1,• • 


(2-40) 


Recalling the nonincreasing property Eq. (2-7) of the diagonal elements 
r^i/ it follows from Eq. (2-40) that 


min 
i=o, • 


l 


Equation (2-41) is an indication that column a® is less important to 
the least-squares approximation than the succeeding columns a3, 
j=o+1,••A. Sufficiently small values of the minimum in Eq. (2-40) 
relative to the other |rjjXi| values suggest that a° might be 
eliminated from A. Similar considerations relating to the numbers 
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Jrjjyj|, j=1,•••,£, suggest the elimination of a (possibly different) 
column aY from A. Taken together. Criteria 1 and II may suggest up to 
three distinct candidates for elimination from A. 

2.2.2.2 The Elimination Sequence 

Having specified the following parameters: 

(1) an actuator candidate collection «V; 

(2) a regulation target collection 33 (e.g., {b L0SX , b L0SY }); 

(3) a least-squares weighting matrix W (notational references 
suppressed); and 

(4) a least-squares minimum tolerance parameter to > 0? 
the elimination sequence can be defined. 


Step 1: Initial Factorization . Set A*°* = A = [jl ] and compute a QR- 
decomposition of A* 0 ). [Because of the assumption Eq. (2-33), the 
R-matrix in the decomposition (cf. Eq. (2-3)) 


has the structure of Eq. (2-4c).] 


(2-42) 


Step 2: Initial Least-Squares Solutions . Denote A* 1 * = A*°*P^ and 
solve the least-squares problems of the form Eqs. (2-34) and (2-35), 
respectively, for A^ 1 ). [It follows from Theorem 2-4 that there are 
infinitely many least-squares solutions (cf. Eq. (2-16)), but all give 
zero for the minimum values Eqs. (2-36) and (2-37) corresponding to 
A<1>.] 

Step 3: Initial Elimination . Truncate A^ 1 ) by retaining only its 
first n columns; denote the result by A^2) ; i.e., 

A 12 ’ = [a*! 01 ; .,.:a E 'l ( " 1 ] i n x n 
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where p-| is the permutation mapping of the indices induced by 

the column permutation matrix of Step 2. [Because of the 

assvmptio'n Eq. (2-33) and vhe structure of the matrix R^), the first n 
columns of A^ 1 ) are linearly independent, Moreover, A^) inherits 
the QR-decomposition 


from Eq. (2-42) in which the R-matrix has the form of Eq. (2-4a)'. There 
are unique solutions x^ 2 ) and y< 2 ) to the least-squares problems of 
the form Eqs. (2-34) and (2-35), respectively, for A^ 2 * (cf. Eq. 

(2-17)), and the corresponding minimum values Eqs. (2-36) and (2-37) are 
again zero.] 


Step 4: Secondary Elimination . Examine the numbers r^ j 

i=1,»»»,n; and associated with A^, x^, and y^. 

Using either Criterion I or II, select a column, say the c-th, of A^ 2 ^ 
for elimination. Truncate A^ 2 ^ by removing the o-th column. (When the 
criteria suggest more chan one column of A^ 2 ^ for elimination , all 
eliminations may be done in parallel. The result is a (generally 
expanding) lattice of alternatives for elimination.] 


Step 5: Secondary Factorization . Denote 

a (31 4 [ a Pi(>> ; ...: a Pi(a)-i : a n j. 

and compute a QR-decomposition of A^). [The columns of A^) are 
linearly independent, having been inherited from a( 2). The QR- 
decomposition 
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has an R-matrix with the structure of Eq. U-*b).] 


Step 6; Secondary Least-Squares Solutions . Denote » &^)p(2) and 

solve the least-squares problems of the form Eqs. (2-34) and (2-35), 
respectively, for a( 4), [It follows from Theorem 2-5 that each of the 
two least-squares problems has a unique solution (cf. Eq. (2-18)), and 
the minimum values corresponding to Eqs. (2-36) and (2-37), respectively, 
are generally_jionzero (cf. Eq. (2-19)).] 

Step 7: Subsequent Eliminations . Repeat Steps 4 through 6, 
successively, applying Step 4 of the current cycle to the results of Step 
6 for the preceding cycle. 

Step 8: Termination . Assume that the preceding eliminations have 
reduced the initial collection jtl to the subcollection id' and that each 
elimination (say to u/") from via the application of Steps 4 through 
6 leads to corresponding minimum values (cf. Eqs. (2-36), (2-37)) such 
that either 


or 


(2-43) 

(2-44) 


where A" = [«*/"]. Terminate the elimination sequence;is a minimal 
selection. 

The actual numerical implementation of the termination tests Eqs. 
(2-43) and (2-44) is modified slightly when very small residual vectors 
in the least-squares solutions are required, as in the applications ;o be 
presented subsequently. The test Eq. (2-43) has the form 
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(2-45) 


where p = (pi,»**,p n ) T is the residual vector b L0SX -A"x" at a 
least-squares solution x" corresponding to A". The actual computation of 
the squares of small numbers in Eq. (2-45) in undesirable. In view of 
the equivalence of all norms on E n , and in particular, from Jensen's - 
inequality (Ref. 2-19]: 


max 
i=»l, • 


kf '.<- («, kf)’ /2 1 H, kil £ ■ 


the test Eq. (2-45) is equivalent to the condition 


s kil > l2 - 46 > 

for appropriately chosen ti« Similar remarks apply to the test Eq. 
(2-44). In addition to being simpler and subject to less computational 
error, Eq. (2-46) also provides information (not provided by Eq. (2-45)) 
as to which of the modal coordinates is associated with the unacceptable 
approximation. 

It should be observed that the elimination sequence has embedded 
within it a considerable degree of design flexibility that goes far 
beyond the choice of the parameters jU, 2B , W, and to* Step 3 
accomplishes a (generally quite large) reduction in the number of 
actuator candidates from y (cf. Eq. (2-32)) to n. The candidates that 
are retained in Step 3 are determined by the column permutation pk) 
generated in Step 1. That permutation in turn is the result of 
successive applications of the Householder Lemma 2-2 that are performed 
in Step 1 in order to achieve the diagonal dominance property Eq. (2-6) 
in the factorization of Eq. (2-42). In general, pk) i s no t unique in 
this regard. The implications of this lack of uniqueness have yet to be 
fully explored. 
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2 . 3 Applications with ACOSS Model No. 2 


In the description of the scientific experiment for active control 
synthesis with generic structural models [Ref. 2-4), four principal 
elements of the synthesis process are identified as variables of the 
experiment: 

(1) Design of the basic structure to be controlled; 

(2) Selection of reduced-order models; 

(3) Selection of active transducers; and 

(4) Determination of the controller feedback structure. 

For the work summarized in the present report, the first and last of 
these elements are fixed in order to focus on the other two, especially 
active transducer selection. 

The first element, basic structural design, is a significant 
variable in relation to the generic structural models because of the many 
modifications to original designs. In particular, the original 
structural design for ACOSS Model No. 2 (Ref. 2-20) has undergone three 
major revisions (Ref. 2-21). The choice of Revision 1 fo: 1 the work 
reported here represents a compromise which retains improvements on the 
original design in the modeling of the attachments for certain rigid- 
body elements, while avoiding the extremes of very high inherent 
flexibility (Revision 3) or very high inherent stiffness (Revision 4). 

The last element, controller feedback structure determination, is 
held fixed for the work reported here in order to avoid obscuring th» 
view of whatever improvements in performance may be attributable to 
judicious selections of reduced-order models and of active transduceis. 
The choice of the disturbance-rejection control feedback structure [Re.‘\ 
2-22; Sec. 3.6), which assumes full knowledge of the statistics of the 
disturbance, provides an upper bound for control system performance. 

The disturbance to be accommodated is assumed to be an aperiodic 
disturbance with a constant power spectral density at frequencies from 0 
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to 5 Hz, applied in a fixed direction at nodes no. 37 ard 46 of the 
structural model. Full details of the disturbance definition and 
relevant assumptions have been given previously [Ref. 2-5, Sec. 2). 

To assist in the selection of reduced-order models, criteria for 
ranking the structural modes are commonly used [e.g., Ref. 2-2]. In the 
present context, ranking of modes is done by comparing open-loop modal 
responses of the regulated variable to the disturbance input,* i.e., the 
root-mean-square 4-RMS) value of the response of the two-input, single¬ 
output system between the disturbance input and. the contribution of each 
mode to the line-of-sight error (cf. Eqs. (2-28) and (2-30)). The 
precise formulation has been given previously [Ref. 2-5, Sec. 5]. The 
resultant ranking of flexible-body modes for the chosen structural model 
(Revision 1) is shown in Table 2-1. An inevitable feature of any such 
ranking procedure is that, in general, a selected group G of 
consecutively-ranked modes is not contiguous in frequency; i.e., there 
exist modes not in G which are interlaced in frequency with modes in G. 
For example, the first eleven modes according to the ranking of Table 2-1 
are interlaced in frequency with seven modes (11, 15, 17-21) of 
considerably lower rank. The interlacing can be seen clearly in the 
frequency listing of Table 2-2. A reduced-order model having this 
property is said to be .interlaced ; otherwise i't is called contiguous . 

As indicated in the preceding section, the objective of the active 
transducer selection algorithm (roughly speaking) is to choose a 
relatively small number of those actuators that are expected to have 
relatively high effectiveness in achieving the desired regulation of 
line-of-sight error. It therefore seems desirable that the selection 
process lead to a rich lattice of alternatives (Step 4) for elimination 
from the initial class of actuator candidates. This qualitative property 
is termed richness . It is certainly desirable that upon termination of 
the elimination sequence, the number of retained actuators is 
substantially smaller than the number of modes to be controlled. 

Denoting by p the number of actuators retained in a minimal selection, 
the level of reduction 
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A = A ( J , t n ) 


n - p 


(2-47) 


is an important quantitative measure of the effectiveness of the 
selection algorithm.' It is obvious that the level of reduction should 
depend on the tolerance parameter to» it also turns out to be.a very 
strong function of the class of actuator candidates chosen. 

Principal observations with respect to the performance of the 
selection algorithm that have emerged from work to date are the 
following: 

(1) Level of reduction increases with augmented richness; 

(2) Richness is augmented when: 

(a) Interlaced reduced-order-design models for which active 
transducer selection is tailored are made contiguous; 

(b) The ciass of actua:or candidates is enlarged; 

(c) Suitable restrictions are imposed upon a sufficiently 
large actuator candidate class; 

(3) Minimal selections which realize the same level of reduction 
are not necessarily unique. 

Four design examples with AC0S3 Model No. 2 are described next. 
These examples demonstrate the observations noted and show in addition 
some of the design flexibility that is available when relatively high 
levels of redaction are achievable. 

2.3.1 Specifications for Design Examples 

Considerations which apply to each of the four examples are the 
following: 

(1) The basic structural design is Revision 1 of ACOSS Model 
No. 2 [Ref. 2-21] with Z = diag(O.OOl) in Eq. (2-28); 

(2) A broadband disturbance [Ref. 2-5, Sec. 2) is applied to 
nodes no. 37 and 46 of the structure, assuming constant 
direction; 
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(3) Selection of- modes for inclusion in reduced-order models for 
controller design is assisted by the disturbance response 
ranking of Table 2-1; 

(4) The regulation target collection consists of the elements 
defined-by Eq. (2-31) (whose dimension is that of the design 
model); 

(5) The least-squares weighting matrix W has the form: 

W = diag(w^): n x n (2-48) 

where w^ is the RMS LOS-error disturbance response 
associated with mode of the reduced-order design 
model as listed in Table 2-1; 

(6) The least-squares minimum tolerance parameter is 

t q « 1<f 6 (2-49) 

The four examples are distinguished by differing specifications on 
the reduced-order design model and actuator candidate class «//. The 
specifications are summarized in Table 2-3 and are discussed briefly 
below. 

Reduced-order design model . In Example A, the reduced-order design model 
consists only of the modes ranked 1 through 11 in Table 2-1; it is 
interlaced. In all the succeeding examples, this design model is made 
contiguous !y incorporating the formerly excluded modes between mode 7 
and mode 24. However, mode 19 -- a solar panel torsion mode -- is 
excluded from the numerical implementation. The associated entry in the 
least-squares weighting matrix W (cf. Eq. (2-48)) is extremely small 
(10“ 22 ) causing the resulting weighting matrix to be poorly conditioned 
numerically with respect to inversion. The small weighting is a 
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reflection of the very low influence of this mode (ranking at most 147 
out of # 156) on any of the modeled node-connecting elements. 

Actuator candidate class ,?l . In Examples A and B, only axial member 
actuators (137 total) on node-connecting elements of the finite-element 
model are allowed. In Example C, the class of candidates is expanded to 
include translation actuators at all nodes, an-aggregate of 311. In 
Example D, certain nodal translation actuators are disallowed; in 
particular, these are the z-translation actuators on equipment section 
nodes (6 total) and both y-translation and z-translation actuators on 
interior solar panel nodes (16 total). 

2.3.2 Selection Results 

Active transducer selection results with Examples A through D are 
summarized in Table 2-3; visual representations are given in Figures 2-7 
through 2-9. The results are discussed briefly. 

Example A yields somewhat unsatisfactory results from the 
viewpoint of design flexibility in actuator selection (A=1)• Only two 
alternatives are available (Step 4) for removing an actuator from the 
candidates remaining after completion of the initial elimination (Step 
3). Both alternatives lead to least-squares solutions which exceed the 
specified termination tolerance (Eq. (2-49)). Hie least undesirable of 
the two elimination alternatives leads to the minimal selection shown in 
Fig. 2-1. It can be seen that the measure of the residual vector in Eq. 
(2-46) slightly exceeds the specified tolerance to in the regulated 
varible associated with b L0SX (cf. Eq. (2-31)). Inclusion of the 
member connecting one solar panel (at node 52) to the equipment section 
(at node 43) is somewhat surprising, since in a static comparison of 
influence upon each of the modes in the design model (intra-rcw 
comparison of entries in the $ T B^ matrix of Eq. (2-28)), this member 
ranks no higher than 91st of all 137 members. 
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Example B yields a somewhat richer lattice of elimination 
alternatives and a substantial improvement of the level of reduction in 
the minimal selection. This is attributable entirely to having 
(minimally) augmented the reduced-order design model to make it 
contiguous, as Table 2-3 indicates. The minimal selection is shown in 
Fig. 2-2. 

Example C yield, only a slightly richer lattice of elimination 
alternatives, but a quite substantial improvement in the level of 
reduction in the minimal selection, in comparison to example B. This is 
attributable to having enlarged the class of actuator candidates to 
include translation actuators at nodes, as Table 2-3 indicates. This 
explains the presence of in Eq. (2-47). The minimal selection is 
shown in Figure 2-3. It is noteworthy, but not surprising, that no axial 
actuators appear in the minimal selection even though the candidate class 
includes them all. In fact, all axial candidates are removed in the 
initial elimination (Step 3). 

Example D provides somewhat surprising and quite satisfying 
results. The lattice of elimination alternatives is extremely rich and 
leads to an even more substantial improvement in the level of reduction 
in minimal selections, of which there are two in this case. A portion of 
the branch of the elimination lattice that leads to the minimal 
selections is shown in Fig. 2-4. The label "TX-34" indicates the removal 
of the translational (T) actuation in the x-coordinate (X) body-axis 
direction at node 34, other labels having an analogous interpretation. 
Parentheses indicate eliminations along a portion of the branch that do 
not lead to one of the minimal selections. Preselections D-17, D-14, and 
D-11 are shown in Figs. 2-5 through 2-7, respectively, and minimal 
selections D-9A and D-9B are shown in Figs. 2-8 and 2-9, respectively. 
Taken together, Figs. 2-4 through 2-9 show quite clearly how the 
elimination sequence (beginning with Step 4) leads systematically to 
minimal selections. The steady increase in the measure of Eq. (2-46) for 
the residual vectors associated with b L0SX and b L0SY (Bq. (2-31)), 
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respectively, as the elimination sequence proceeds may be noted, 
consistent with the expectation of Theorem 2-6 (Eq. (2-27)). In this 
example, the number of actuators in a minimal selection is only half the 
number of controlled modes. This implies a very desirable degree of 
design flexibility is available (and may be beneficial) by employing 
intermediate (ancestral) selections such as D-14 or D-11 that lie along 
the p'ath to a minimal selection.' A full explanation of how the 
restriction imposed in proceeding to example D leads to the improvements 
observed is net yet available. However it does provide a substantiation 
of the design flexibility embedded within the elimination sequence that 
is mentioned at the end of Section 2.2. 

Minimal selections from Examples B, C, and L are employed in the 
design oJ: optimal linear-quadratic disturbance-rejection controllers as 
reported subsequently in Section 3. Each of the designs is evaluated in 
the presence of a family of contiguous reduced-order evaluation models 
(cf. Table 2-2) with increasing ranges of frequency.- Stability well 
beyond the frequency range (5 Hz) of the design model is demonstrated in 
each case. 

2.4 Summary 

Work accomplished to date in connection with the scientific 
experiment for active control synthesis with broadband disturbances has 
led to substantial insight into a relatively new and systematic approach 
to active transducer selection. This process leads to selections which 
at the same time have relatively high influence upon variables to be 
regulated in the control problem and require substantially fewer 
actuators than controlled modes. The selection process itself, as well 
as the selections it produces, exhibit potentially beneficial design 
flexibility. Linear-quadratic disturbance-rejection controllers 
employing transducer selections generated by this process (cf. Sec. 3) 
demonstrate stability in the presence of residual modes over a frequency 
range substantially beyond that of the design modes. These 
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demonstrations are significant in that stability-enhancing adjustments of 
the controller feedback structure have been deliberately postponed in 
order to focus on the connections between the. selection of reduced-order 
models and active transducers on the one hand and closed-loop stability 
and performance on the other hand. 

Subsequent work in regard to active transducer selection will 
include: 

(1) Expansion of the actuator candidate class to include 
rotational actuators; 

(2) Deeper investigation of the design flexibility inherent in 
the selection process; and 

(3) Examination of the interaction between the selection of 
reduced-order models and the selection of active transducers. 

2.5 Appendix 

Theorem 2-6 is an immediate consequence of the following two 
results. 

Proposition 2-7 (Column Rearrangements).. The minimum value defined by 
Eq. (2-14) is unchanged by an arbitrary column permutation P: mxm of Aj 
i.e., 

"^ n J(P _1 x ; AP, b) = II ^ n J(x; A, b) (2-50) 


Proof . The identity 

J(P _1 X; AP, b) = Jb-(AP)(P -1 X)J 2 = J(x; A, b) (2-51) 
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holds for all n-vectors x. It follows that for each minimum, say x 0 
(cf. Eg. (2-13)), of J(•; A, b), P^Xq is a minimum of J(•; AP, b), 
and conversely. The desired coincidence of minimum values Eq. (2-50) 
follows from Eq. (2-51) by this correspondence of minimum points. □ 

Proposition 2-8 (Terminal Removals). Let p be an integer, 1 £ p £ m-1. 
The minimum value defined by Eq. (2-14) does r.ot decrease if columns 
p+1,»»»,m are removed from A; i.e., 


m J(xj A, b) < 


min ... 

„p J(y? A. ; 
yeE F 11 1...p' 


where Aidenotes the submatrix consisting of the first p columns 
of A. 


Proof . Integers k, m, and n are defined as in Eq. (2-1). In view of 
Proposition 2-7, there is no loss of generality in assuming.that the 
leading k columns are linearly independent. Two cases are considered, 
depending on whether Eq. (2-1) represents an equality or a strict 
inequality. 

Case A; k = min{n,m> . If n £ m, then the result is trivial: Eq. (2-15) 
is true, and therefore 


J . a 0 < J P . 

nun — mxn 

Thus, consider that m < n. The QR-decomposition of A (cf. Theorem 2-1) 
has the structure 


A - ® 5 : Q 2 ] 


(2-52) 
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where in this case the column permutation P is taken to be the identity 
matrix. The truncated matrix Ai...p inherits a QR-decomposition from 
Eq. (2-52) by a corresponding truncation of R: 


V..p - ®1...P s fc:«2] 

where the matrix Q is the same as in Eq. (2-52) but is repartitioned to 
be compatible with (Ri)i... p : pxp. Denoting by q3 the j-th column 
of Q, it can be seen that 


Using Eq. (2-19) with both $ 2 and Q^, it follows that: 


jP 

min 


i j=P+1 




'«j b ' ^ J »i 


Case B: k < min {n ,m} . The proof for this case uses the same approach 
as in Case A, except that Eq. (2-52) is replaced by the decomposition for 
the rank-deficient case that is a simple extension of Theorem 2-1: 


A = QR 5 [Q, : Q 2 ] 



The details are omitted. □ 
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Table 2-1. Ranking of modes: disturbance input to open-loop modal 
RMS LOS-error output; ACOSS Model No. 2, Revision 1. 
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Table 2-2. Reduced-order models in terms of open-loop 
undamped characteristic frequencies: 

o Interlaced design 
• Contiguous design 
■ Contiguous evaluation 


MODE EIGENVALUE REAL EIGENVALUES 

NO. J2 <*>l w, 

1 (rad/s) (Hz) 


1 

2 

3 

4 

5 

6 




0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 

0 

• 

■ 

8.650672E-01 

9.300899E-01 

1.480284E-01 

8 

O 

• 

■ 

3.148905E+00 

1.774515E400 

2.824226E-01 

9 

0 

• 

■ 

4.009007E+00 

2.002251E4Q0 

3.186681E-01 

10 

o 

• 

■ 

4.434997E400 

2.105943E400 

3.351712E-01 

11 


• 

■ 

8.638482E*00 

2.939130E400 

4.677770E-01 

12 

0 

• 

■ 

1.344033E+01 

3.666106E400 

5.634789E-01 

13 

0 

• 

■ 

1.424241E401 

3.773911E400 

6.006367E-01 

14 

0 

• 

■ 

1.79032&E+01 

4.231224E400 

6.734203E-01 

15 


• 

■ 

3.639226E401 

6.032599E40C 

9.601181E-01 

16 

0 

• 

■ 

4.708934E+01 

6.862167E400 

1.092148E400 

17 


• 

■ 

1.334769E*02 

1.1S5322E401 

1.8387S2E400 

18 


• 

■ 

1.341718E+02 

1.158326E401 

1.843533E400 

19 


• 

■ 

1.40895SE*02 

1.186994E401 

1.889159E400 

20 


• 

■ 

1.564057E402 

1.250623E401 

1.990426E400 

21 


• 

■ 

1.675513E+02 

1.294416E401 

2.060127E400 

22 

O 

• 

■ 

2.372882E+02 

1.540416E401 

2.451648E400 

23 

O 

• 

■ 

2.412209E+02 

1.553129E401 

2.471881E400 

24 

O 

• 

■ 

4.149673E*02 

2.037074E401 

3.242105E400 

25 



■ 

1.051814E403 

3.243167E401 

5.161660E400 

26 



■ 

1.055404E403 

3.248690E401 

5.170464E400 

27 



■ 

2.449482E403 

4.949223E401 

7.87&933E400 

28 



■ 

2.474952E403 

4.974889E+01 

7.917701E400 

29 



■ 

3.037863E403 

5.511682E401 

6.772115E400 

30 



■ 

3.040462E403 

5.514038E+01 

8.775864E400 

31 



H 

3.083413E403 

5.552849E401 

8.637633E400 

32 



■ 

3.183590E403 

5.642339E401 

8.9D0061E400 

33 



■ 

3.653264E403 

6.044223E401 

9.6196796*00 

34 



S 

4.265980E403 

6.531447E401. 

1.039512E*01 

35 



■ 

5.379633E403 

7.334598E401 

1.167338E»0L 

36 



■ 

5.939781E403 

7.706998E401 

1.226607E401 

37 



■ 

7.075340E403 

• 8.411504E401 

1.338732E*01 

38 



■ 

7.319027E403 

8.S55132E401 

1.361592E*01 

39 



■ 

8.614004E403 

9.281166E401 

1.477143E+01 

40 




1.065900E+04 

1.032424E+02 

1.643153E401 

41 




1.713255E404 

1.308914E+02 

2.083200E401 

42 




1.674883E404 

1.369264E402 

2.179250E401 

43 




1.883B82E404 

1.372546E402 

2.134474E401 

44 




1.884647E404 

1.372824E402 

2.184917E*01 

45 




1.884669E404 

1.372832E402 

2.184930E*01 

46 




1.889893E404 

1.374734E402 

2.187956E*01 

47 




2.227858E404 

1.492601E402 

2.375548E401 

48 




2.364327E404 

1.537637E402 

2.447224E*01 

49 




2.455703E404 

1.567068E402 

2.494066E401 

5(1 




2.643350E404 

1.625838E402 

2.587601E401 
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Active transducer selection: design example specifications and results. 













































SECTION 3 


PROGRESS ON SYNTHESIS OF ACTIVE CONTROL FOR BROADBAND 
DISTURBANCE ACCOMMODATION 

PART 2s CONTROL SYSTEM DESIGN AND EVALUATION 


3.1 Introduction 

The philosophy of our scientific experiment for active control 
design with ACOSS Model No. 2 was reviewed in the previous section. Now 
we make use of the 18 mode (modes 7-24) design model from Revision 1 of 
Model No. 2, and let the several actuator and sensor placements be 
variables in the experiment. 

This section will review the control system design philosophy 
which remains fixed while working with the various actuator and sensor 
selections. After the design results are examined, some changes in the 
controller design will be suggested for the future part of the 
experiment. 

3.2 Control System Design Philosophy 

For these designs we have employed the Linear-Quadratic-Gaussian 
(LQG) methodology, where a model of the disturbance has been augmented 
with the state equations of the structure. Such a disturbance model uses 
knowledge of the point of application on the structure and power spectral 
density of the incoming disturbances. The resulting controller is 
frequently referred to as a disturbance rejection controller, and has 
been treated in References 3-1 (Section 4) and 3-2. 
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To clarify some of the values used here and in the designs of 
Reference 3-1, we examine the transfer function model for the power 
spectral density function corresponding to the incoming disturbances as 
given by Eq. (2-1) of Reference 3-1 (Section 2): 


where 



G 


40 N 2 
2irf 


c 


sec 


5 Hz 


The corresponding state space representation is (cf. Eq. (2-28)) 

V = BqUd 

where uq is the solution of 


u_ = JLu + B 0 ) 

D u D ID 

In particular, v is the disturbance acting on the structure, and b d , 
the influence matrix of the disturbances, has been discussed in Section 2 
of Reference 3-1. Here 


where Ai = X 2 ~ -uic = -31.4159, and 
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ail 


.0 


b 


0)2-1 


where b^-j = b^ = /<3 uic = 198.692. The above terms are discussed 
in Section 4 of Reference 3-1, but without noting the numerical values. 
Also the covariance matrix of the fictitious white noise input vector m 
is 


Cov(oi ai T ) = 12x2 

Note that the power spectral density function has a bandwidth of 5 
Hz and has a zero frequency value of 40 N 2 /Hz if it is considered to be 
a two-sided function of frequency, and has a zero frequency value of 80 
N 2 /Hz if it is considered to be a one-sided function of frequency. 

Among the controller designs presented in this section for the 
four actuator/sensor placements, the only parameters that have been 
altered are the diagonal terms of the 3x3 Q matrix in the performance 
index for the LQ regulator where, 

’q, 0 0 ' 

Q * 0 q 2 ° 

■° 0 * 3 . 


and q-j weights the line-of-sight displacements in the x direction 
(LOSX), q 2 weights the line-of-sight displacements in the y direction 
(LOSY), and q 2 weights the defocus displacements (Z). In particular, 
the diagonal terms of the Q matrix are adjusted to ensure that the stare 
feedback portion of each LQG design will meet the line-of-sight (LOS) 
specifications. These LOS specifications are taken to be: 
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® 1 /2 

LOSX = [ £ (Oi - Y1 )] ' <_ 0.05 jj- radians 

i=’ 

^ - /2 

LOSY = [ l (Ui - Y2)] ' £ 0.05 y-:radians 

i=1 

Defocus = Z = [ £ (Ui - Y3)] 1 / 2 <_ 0.025 x 10 3 meters 

i=1 

where Ui-Yj represents the mean square «alue at one of the three line- 
of-sight outputs resulting from the i-th disturbance. For our designs 
there are two identical disturbances, one at node 37 and one at node 46 
of the model, so m is two in the above LOS specifications. 

For all designs the control weighting matrix, R, in the LQ 
performance’ index was an identity matrix of appropriate dimension and its 
scalar factor, p R , was always 1x10“ 17 . 

The design of the optimal state estimator is dual to that of the 
LQ regulator and the covariance of the sensor noise term, 9, was taken to 
be: 


Covfe 0 T ) = pi 

' y o mxm 

where p Q = 1x10 and m is the number of sensor outputs, as in the 
designs of Section 4 of Reference 3-1. 

The state estimator is connected to the constant LQ gain to form a 
compensator which, when evaluated with the 18 mode design model of the 
structure, is guaranteed to yield a stable closed-loop system. Note that 
Model No. 2 has 132 elastic modes of higher frequency than the 18 modes 
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used in our design. One goal should be to have a control system robust 
enough to remain stable in the presence of these 132 high frequency 
modes, realizing that degradation in LOS performance would occur. Our 
designs have not achieved that goal at this point in the experiment, but 
the results, discussed in the next section, are encouraging. 

3.3 Controller Results 

The results of the control system designs with four of the 
actuator/sensor selections described in Section 2 are summarized in Table 
3-1. The diagonal elements of the Q matrices used in the LQ designs, and 
the resulting line-of-sight results for the full-state feedback case 
under the column headed by LQ are shown. Note that the LQ designs with 
translational actuators/sensors required much lower weighting of the 
controlled variables (line-of-sight displacements) in the performance 
index than did the design with axial members, while achieving 
approximately the same LOS. performance. 

Not reflected in the table are the design iterations on the Q 
matrix for each particular actuator/sensor selection that led to 
acceptable LOS performance. These iterations showed that each Z result 
was relatively insensitive to changes in its own weighting or the 
weighting of the other terms. However, an order of magnitude change in 
the weighting for LOSX or LOSY, namely qi or q 2 , from those shown in 
the table yielded approximately an order of magnitude change in the LOSX 
or LOSY value. 

Recalling that actuator/sensor selection D-9A was a subset of 
selection D-14, one can note only a slight degradation in the LOS 
responses between these two for the LQ design. 

The original 18 mode LQ design for selection D-9A encountered 
numerical difficulties when evaluating the poles of its state-feedback 
design. These numerical difficulties were alleviated when mode 19, which 
has extremely low control and output influences, was deleted from the 
design process. Mode 19 was of course included in the evaluation models. 
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Going to the final column of Table 3-1 , one can see that 
actuator/sensor selection C yielded the largest number of modes in a 
stable evaluation. The addition of one mode beyond any of the stable 
evaluations shown in the last column gave an unstable system. Note that 
selection D-9A, which was a subset of selection D-14, had a much lower 
stable evaluation model than D-14. In general, taking into account the 
number of actuators/sensors in each selection, the translational members 
yielded control systems with superior stability properties compared to 
the design using axial members. 

For the two designs giving the greatest envelope of stability, the 
LOS results were computed for the LQG compensator evaluated with the 18 
mode design model and with a 30 mode evaluation model (modes 7-36). 

While slight degradation occurred in the LCSX and LOSY results for the 
LQG compensator evaluated with 18 modes compared to the full-state 
feedback results in the LQ column, significant degradation in LOS did 
occur when the LQG compensator was evaluated with modes 7-36. 

3.4 Future Controller Designs 

As part of the scientific experiment explained in Section 2, 
changes will be made in the control system design with an eye toward 
improving the stability results. Specifically, the choice of state and 
control weighting matrices in the LQ performance index could be altered 
to lower the multivariable bandwidth. In addition, other frequency 
shaping techniques may have to be used to ensure adequate high frequency 
attenuation in the LQ design. Hopefully, in addition, a reasonable 
disturbance rejection controller design can be found which relies less on 
explicit knowledge of the disturbances' structural locations and 
stochastic properties. If a linear quadratic regulator can be designed 
with some of the improved properties mentioned above, the estimator may 
have to be designed using "full-state recovery techniques" [Ref. 3-3) to 
preserve those good features, which otherwise are not guaranteed to hold 
when an estimator of general design is used in the compensator. 
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1.5 


Conclusion 


This section has illustrated the LQG disturbance rejection 
controller results for four actuator and sensor selections used on ACOSS 
Model No. 2. These selections were the variables under investigation in 
the experiment referred to in Section 2. The results indicated that 
translational actuators and sensors gave superior line-of-sight and 
stability results when compared to axial actuators and sensors. Changes 
in the controller design philosophy will be made to improve the stability 
results in the future part of this experiment. 
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SECTION 4 


SAMPLED DATA CONTROL OF FLEXIBLE STRUCTURES 
USING NON-COLOCATED VELOCITY FEEDBACK 


4.1 Overview 

A framework is developed for sampled data control of flexible 
structures, in terms of discrete time recursive equations in second order 
form. This framework is used to analyze the sampled data control scheme 
where the loop is closed using constant gain output velocity feedback. 

It is well known that the closed loop is stable if colocated velocity 
feedback with symmetric and positive definite feedback gain is used, so 
long as the sampling rate is sufficiently high. In this section it is 
shown that the closed loop can be stabilized using sampled data output 
velocity feedback for arbitrary sampling rate; our approach is to make 
use of appropriately defined non-colocated velocity feedback. This 
approach leads to explicit stability conditions in terms of the feedback 
gain and the sampling rate. 

4.2 Models for Sampled Data Controlled Flexible Structures 

A sampled data controlled flexible structure can be defined as a 
distributed parameter system, where the structure input is the output of 
an ideal zero order hold and the structure output is sampled. Although 
distributed parameter models typically involve infinite dimensional 
variables, our analysis is based on the finite dimensional model 

Mq + Kq = Bu (4-1 ) 
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For simplicity in the subsequent development no structural damping is 
included. The structural displacement vector q = (q-j, • •• ,q n ) and the 
force input vector u = (ui, • • *, 1 ^). The mass matrix M and the 
structural stiffness matrix K are assumed symmetric and positive 
definite. Throughout, we consider velocity output of the form 

y = Cq (4-2) 

where the velocity output vector y = (y-j, • • •,y m ). The input 
influence matrix B and output influence matrix C are assumed to be 
dimensionless. 

The structure input u is defined in terms of the input sequence 
u^ by the ideal zero order hold relation 

u(t) = u^ , kT <_ t < kT+T (4-3) 

The output sequence y* is defined in terms of the structure, output y by 
the ideal sampling relation 


y k = y(kT) 

The fixed value T > 0 is the constant sampling time. This open loop 
sampled data controlled structure can be viewed as a discrete time system 
with input sequence u^ and output sequence y^, where k = 0,1,• 

Let $ be a nxn nonsingular modal matrix and let f) 2 be a nxn 
diagonal modal frequency matrix (Ref. 4-1] satisfying 

<t T M$ = I, 4>' f K$ ■ A 2 = diag(oy •••,<!)*) 


Introduce the coordinate change 


4-2 





q = $n 

so that Eqs. (4-1), (4-2) can be written as 

n + fl 2 n = * t bu (4-4) 

y = ct n 

It is an easy task to solve the vector equation Eq. (4-4), using Eq, 
(4-3), to obtain 


\+1 = (cosOT) \ + n” 1 {sinflT)^ + n" 2 (i - cos£)T) (4-5) 

n k+1 = -nfsinflT)^ + (cosQTjf^ + fl' 1 (sinftT) (4-6) 

where n k = n(kT), = n(kT), and sinflT = diag( sin( T), • •*,sin( u^T)), 
cosfiT = diag(cos(<uiT), • • ♦,cos((u n T)) (Ref. 4-2]. 

Although the first order recursive equations Eqs. (4-5), (4-6) 
could be used, it is more convenient for our purposes to make use of a 
second order recursive equation in n k alone 

n k+1 - 2(cosnT)n k + n k-1 = iT*(sinnT)$ T B(u k - u Jt _ 1 ) (4-7) 

y k = C$ n k (4-8) 
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The modal equations Eqs. (4-,7), (4-8) form the basis for our subsequent 
analysis. It is natural to make use of the recursive equation for 
alone in considering velocity feedback systems. 

It should be noted that relations Eqs. (4-7), (4-8) involve no 
rtumerical approximation; they are valid for any sampling time T > 0. 

4.3 Constant Gain Velocity Feedback 

Constant gain output velocity feedback has been studied 
extensively for analog controlled structures. Our interest is in use of 
constant gain output velocity feedback for sampled data controlled 
structures. 

Consider the closed loop sampled data controlled structure defined 
by Eqs. (4-7) and (4-8), using the control input sequence 

\ - -Gy k (4-9) 

where G is a constant mxm feedback gain matrix. Substituting Eq. (4-9) 
into /Eqs. (4-7) and(4-8), a closed loop recursive equation is obtained 

n k+1 - [2cosOT - fl" 1 (sin«T)<S> T BGC<&]n k 

+ [i - n' 1 (sinam T BGc*]n k-1 - o (4-io) 


The closed loop characteristic eq ■ • .on can be written as 

= 1 - (2cosftT - ft” 1 (sinftT)$ T BGC<&)z 

/<I - ft” 1 (sinftT) $ T BGC$) ] = 0 


d(T,z) 









The objective of constant gain velocity feedback control is to make the 
closed loop as described by Eq. (4-10) geometrically stable, i.e., to 
make the closed loop characteristic zeros lie inside the unit disk. 

We use Eq. (4-10) as the basis for our subsequent analysis of the 
closed loop. If sinflT is nonsingular, the following implications hold: 
if n^+0 as k+«, then necessarily u^O as k+«> and n^+0 as k +<*>} 

consequently x^+0 and x^+0 as k+«. 


4.4 Colocated Velocity Feedback 

Recall the following results for constant gain output velocity 
analog feedback control where u ■ -Gy. If colocated force actuators and 
velocity sensors are selected so that C = B T , then the closed loop 
analog controlled structure is asymptotically stable if G is any 
symmetric, positive definite matrix, and if a certain controllability 
assumption is satisfied [Refs. 4-3, 4-4, 4-5]. Moreover, this result 
does not depend on the particular values of the modal frequencies and 
modal functions. 

We first mention a rather obvious result that if the sampled d?ta 
feedback control is chosen according to the analog feedback theory, then 
the closed loop is stable for sufficiently small sampling time. The 
brief proof is included for completeness; it also serves as an 
introduction to our subsequent development. 

T heorem 4-1 . Assume that 
. (a) C = B t , 

(b) the matrix pair ((l 2 , $ T B) is completely controllable, and 

(c) G is symmetric and positive definite. 

Then the closed loop equation Eq. (4-10) is geometrically stable for 
sufficiently small sampling time T > 0. 





Proof . Consider the associated polynomial 


p(T,w) = det [2(1 + cosflT - ff 1 (sinfiT)* T BGC$)^p- 

+ 2JT 1 (sinftT) <& T BGC$ + 2(1 - cosflT)] (4-11) 

3y assumption the polynomial defined by 


lim p(T,w) — 

T+0 T 2 

has all zeros in the left-half-plane; hence there'is a number T > 0 such 
that p(T,w) ha'’ all its zeros in the left-half-plane for 0 < T < T. 

Using the bilinear transformation 


1 + — 

z . ■ (4-, 2 ; 

" 2 

it follows that the zeros of d(T,z) are necessarily inside the disk; 
hence Eq. (4-10) is stable. □ 

This result has limited application since there is no indication 
of the range of values of the sampling times, relative to the feedback 
gain matrix, required for closed loop stability. In Reference 4-6, 
conditions are developed which, in principle, characterize a range of 
values of the sampling time for which the closed loop system is stable. 
Unfortunately, the conditions depend on an a priori computable bound on 
an exponential matrix; computation of such a bound, in analytical terms, 
is not considered in Reference 4-6. Of course one could perform a 
numerical search, based on the characteristic polynomial d(T,z) (or 
equivalently p(T,w)) for a specific case, to determine a range of values 




of the sampling time for which the closed loop is stable. However, for 
the case of colocated velocity feedback there are no known explicit 
conditions, in terms of the sampling time and feedback gain matrix, which 
guarantee stability of the closed loop sampled data system. 

4.5 Non-colocated Velocity Feedback 

We now present the main result of the papers a set of explicit 
conditions on the sampling time and feedback gain matrix for which the 
closed loop sampled data controlled structure is stable. The key idea is 
to suitably modify the assumption of colocated force actuators and 
velocity sensors. 

Theorem 4-2 . Assume that sinftT is nonsingular and 

(a) C = B T <MsinflT)n' 1 T" 1 <r 1 , 

(b) the matrix pair .([I + cosQTj^H - cosfJT], <fr T B) is 
completely controllable, 

(c) G is symmetric and positive definite, and 

(d) I + cosflT - T* t C t GC<S is positive definite. 

Then the closed loop equation Eq. (4-10) is geometrically stable. 

Proof . The assumptions can be shown to guarantee that the zeros of 
p(T,w) defined in Eq. (4-11) are in the left-half-plane. The bilinear 
transformation defined in Eq. (4-12) guarantees that the zeros of d(T,z) 
are necessarily inside the unit disk. Hence Eq. (4-10) is stable. □ 

Informally, note that as the sampling time satisfies T+0 that from 
condition (a), C+B T , and that conditions (b), (c), (d) of Theorem 4-2 
tend toward conditions (b), (c) of Theorem 4-1. But for T > 0 the extra 
requirement of condition (d) of Theorem 4-2 is imposed on the feedback 
gain; this imposes a limit on the "size" of the gain matrix. Condition 
(a) of Theorem 4-2 implies that the output velocity feedback given by 
Eq. (4-2) be selected in a specific way. Clearly the actuators and 
sensors would generally not be colocated. 
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Conditions (a) and (d) of Theorem 4-2 do depend on explicit 
knowledge of the modal frequencies and modal functions and they depend on 
the particular value of the sampling time. Unlike the corresponding 
analog case, the closed loop, with assumptions of Theorem 4-2 satisfied, 
is not robust to arbitrary uncertainties in the modal data and to modal 
truncation. However, an arbitrary structure can be stabilized using 
sampled data velocity feedback with arbitrary sampling time; the 
suggested approach depends on a proper choice of force actuator and 
velocity sensor locations so that condition (a) of Theorem 4-2 is 
satisfied. 

4.6 An Example • 

As a simple example consider the scalar equations in modal 
coordinate form with one actuator input 

"i + Vi - “ 

- 2 

n 2 + 0> 2 n 2 = u 


2 2 

where 

Based on Theorem 4-1, the sampled data feedback control 

u(t) «* -Gfn^kT) + h 2 (kT)] , kT <_ t < kT+T, 

with G > 0 and u^T * nir, n=0,!,•••; i=1,2 stabilizes the closed loop 
for T > 0 sufficiently small. Even in this simple case analytical 
conditions on the range of values of the sampling time and the feedback 
gain for wIJLch the closed loop system is stable are exceedingly 
complicated. Evidently a proper choice of values for the feedback gain 
and the sampling time would require use of numerical search procedures. 
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Now consider the sampled data control 


sinakT sinui_T 

u(t) = ^~^r ~) V kT) + ) V kT) ]' w < t<kT+T ' 

where the output is chosen so that condition (a) of Theorem 4-2 is 
satisfied; clearly such velocity sensors would not be colocated with the 
specified force actuators. On the basis of Theorem 4-2-the closed loop 
is guaranteed to be stable if o^T * nir, n=0,1,«»*; i=1,2; and if the 
relatively simple analytical conditions 

(a».T) 2 (1 + costu T) 

0 < GT < —---—-— , i=1,2 

(sinWjTj 

0 < GT < {a» 1 T> 2 (o) 2 T) 2 (1 + cosaijTHl + coso^T) 

(oj 2 T) 2 (1 + cos« 2 T) (sin^T) 2 + (w 1 T) 2 (1 + cosu^T) (sina^T) 2 

are satisfied. Thus proper choice of values for the feedback gain and 
the sampling time is considerably simplified. 

4.7 Conclusions 

We have presented one approach where feedback gains and sampling 
time can be readily selected to guarantee that the sampled data 
controlled structure is stable. Our conditions imply that force 
actuators and velocity sensors be located to satisfy a specific 
analytical condition which, in fact, depends on the sampling time. Ihese 
explicit results are in contrast to the lack of available explicit 
guidelines for choosing the feedback gains and a sampling time in the 
case where colocated actuators and sensors are used. 
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SECTION 5 


OPTIMAL TRACKING AND TERMINAL TRACKING MANEUVERS 
FOR FLEXIBLE SPACECRAFT 


5.1 Introduction 

The problem of feedback control of flexible spacecraft undergoing 
large-angle maneuvers is a topic of continuing interest in the aerospace 
community. Two types of closed-loop slewing maneuvers are of particular 
interest. First, tracking maneuvers where the control is determined in 
such a way as to cause the state to track or follow a desired output 
state. Second, terminal tracking maneuvers where some elements of the 
state are required to exactly satisfy nonzero terminal constraints. 

The necessary conditions defining the solutions for tracking and 
terminal tracking problems are well known. Nevertheless, the 
computational burden required to numerically integrate the resulting 
systems of coupled nonlinear Riccati-like differential equations has 
seriously hindered the common application of these control techniques. 

To overcome the computational difficulties noted above, closed 
form solutions are presented for the necessary conditions defining the 
solutions for the following linear time-invariant control problems: 

(a) The linear tracking problem? 

(b) The terminal tracking problem; 

(c) The state trajectories for the conventional closed-loop 
system dynamics equation. 

The remainder of the section is presented in three parts. Section 

5.2 presents closed form solutions for the Riccati and prefilter 
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equations of the linear tracking problem. Maneuver simulations are 
presented for both conventional and control-rate-penalty control designs 
[Refs. 5-1, 5-2]. Section 5.3 presents closed form solutions for the 
three nonlinear coupled Riccati-like differential equations, defining the 
solution for the terminal control problem. Maneuver simulations are 
presented for both conventional and control-rate-penalty control 
designs-. Section 5.4 presents a closed form solution for the state 
trajectories of a conventional feedback control system, when the plant 
and the state are assumed to be perfectly known. The solution for the 
state is presented in a simple recursive form. 

5.2 Closed Form Solution for the Linear Tracking Problem 

In this subsection the problem of maneuvering a flexible 
spacecraft through a large angle is considered, where the feedback 
control system is required to track a desired output state. The desired 
output state is assumed to be provided from an open-loop solution for the 
linear time invariant system model (see References 5-1 and 5-2 for 
examples of suitable open-loop solutions). 

The optimal control problem of this section is formulated by 
defining a performance index which consists of an integral of quadratic 
forms in the state, control, and control-rate. The necessary conditions 
for this problem lead to two nonlinear Riccati-like differential 
equations [Ref. 5-6]. One equation is the standard Riccati equation, 
which possesses a well-known solution in terms of a steady-state plus 
transient term (Ref. 5-1]. The remaining vector equation is a prefilter 
equation which has foreknowledge of the desired output state, and is 
coupled to the standard Riccati equation. In Section 5.2.4, a new closed 
form solution is presented for the auxiliary Riccati-like vector 
equation. 
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5.2.1 Optimal Linear Tracking Problem 

The optimal reference tracking problem is formulated by finding 
the control inputs u(t) to minimize 

t f 

•J = ‘1 nz*(t f ) - z(t f )ig+-j J [iz*(t) - z(t) + iu(t) n^]dt 

fc o ' (5-1) 


for the system 


x = Ax + Bu , x(t Q ) = x 0 (5-2) 

z = Cx (5-3) 

T 

where lpl w = p Wp, z* is the desired output state, A is the system 
dynamics matrix, B is the control influence matrix, C is tne measurement 
influence matrix, S is the terminal state weight matrix, Q is the state 
weight matrix, R is the control weight matrix, x is the state, and u is 
the control. (For maneuvers where the state is augmented by the control 
and control-rate penalties, the A and B matrices of Eq. (5-2) are 
modified as shown in Reference 5-1). 

5.2.1.1 The Open-Loop Control Problem for the Desired Output State 

The desired output state, z*(t), in Eq. (5-1) is assumed to be 
provided by the solution for the following open-loop control problem, 
where we seek the control inputs u*(t) to minimize 

J* = j J [lx*(t)|^ + Ju*(t)I^]dt (5-4) 

fc o 
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for the system 


x* = Ax* + Bu*, x*(t Q ) 
z* = Cx* 


As shown in Reference 5-2, the closed form solution for the 
open-loop problem defined by Eqs. (5-4) and (5-5) can be written as 


r 

.. \ 

X* 

° \ 

U*(t)J 

1 

hi 


L -Q* 


yv*? - ♦ x oi' 

\ x*A* / l ,x*x* ) 


= the exponential matrix, 


_A*(t) = the costate vector for the open-loop 
solution, 


and <(>_ 


x*x * 1 Vl* 


are partitions of e^ *) for t=t^. 


To obtain x*(t) from Eq. (5-6), the following mapping equation is 
defined: 
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x*(t) 


(5-9) 



where H = [i o] is a (n x 2n) selection operator and x* is (n x i). 

As shown in what follows, the explicit presence of e<•) in Eq. 

(5-9) permits a closed form solution to be obtained for the prefilter 
equation of the linear tracking problem. 

5.2.2 Necessary Conditions for the Linear Tracking Problem 

The necessary conditions defining the optimal tracking problem for 
Eqs. (5-1), (5-2), and (5-3) are gi/en by the following Riccati-like 
differential equations [Ref. (5-3), pp. 100-102 ]: 

•S = -PA - A T ? + PBR -1 B T P - C T QC > P(t f ) = C T SC (5-10) 

l =» -[A - BR -1 B T p] T r - C T QCx* ; 5(t f ) = C T SCx*(t f ) (5-11) 

where the optimal control is given by 

u = -r“'b T [px - ?] (5-12) 

5.2.3 Closed Form Solution for the Time-Varying Riccati Equation 
The solution for Eq. (^-10) has been previously obtained in 

References 5-2 and 5-4, and can be written as 

P(t) = P + 0’\t) (5-13) 

ss 
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where P ss is the solution for the algebraic Riccati equation 


-P A - A T P + P BR~ 1 B T P - C T OC = 0 (5-14) 

ss ss ss ss 

and the closed form solution for 0(t) can be shown to be 

e(t) = 0 ss + a^ (t ' t f ) [0 f - 0 sg ]e A (t "V (5-15) 

where A = A - BR ^^P^, 0 f -(C T SC-P gg ) 1 > an<3 Q gg satisfies the algebraic 
Lyapunov equation 

A0 gg + 0 gg A T => .BR -1 B T (5-16) 

5.2.4 Closed Form Solution for the Prefilter Equation 

The solution for K in Eq. (5-11), follows on assuming the product 
form solution 


Kit) => 0"\t) r(t) (5-17) 

where 0 is defined by Eq. (5-15) and the linear constant coefficient 
vector differential equation for r can be shown to be 

r - Ar = -0 C T QCx* , r(t f ) = 0 f C T SCx*(t f ) (5-18) 

(To derive Eq. >.5-18) one requires the differential equation for 0 given 
by: 0 = A0 + 0A T - BR -1 B T ). 

The solution for «.* in Eq. (5-1 8) can De shown to be 

r(t) e At [r Q - / e _AT 0 ( t )C T QCx*( rl d T ] (5-19) 
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where the initial condition for the vector r is given by 


r Q = e f e f C T SCx*(t f ) + / f e" AT 0(T)C T QCx*(T) d T (5-20) 

o 

The solution for the integral expression in Eq. (5-19) follows on 
introducing 0(t) from Eq. (5-15) and x*(t) from Eq. (5-9) into Eq. 

(5-19), leading to 


/ r AT 0 (t)C T QCx*(t) d T = / e"* x D e® 1 d T 

o L° 



As shown in Reference 5-5, 
exponentials in Eq. (5-21) can be 
constant matrices 



the integrals involving the matrix 
easily evaluated by defining the 
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and computing the matrix exponentials 


’ F 1 (t) 

G^t)* 

V 

‘ F (t) 

G 2 (t) " 



/ e = 



0 

F 2 (t)_ 


0 



where 

F^t) = e St , F 2 (t) = F 4 (t) = e ftt , F 3 <t)- = e - ^ 

«,(« - g 2 (« - 

o o 

As a result, it follows that the integrals in Eq. (5-21) are given by 

/ e’^D^dx = F*(t)G <t> (5-25) 

t -T 

/ e T D 3 e“ T dx = F^(t)G 2 (t) (5-26) 

5.2.4.1 Recursion Relationsnips for the Integrals of Exponential 
Matrices 

In order to evaluate Eqs. (5-25) and (5-26) at discrete time 
steps, the semi-group properties of exponential matrices are now 
exploited; yielding the following recursion relationships for the matrix 
partitions of n C ’ and e C 2<“- 1 >«, 
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F. [(n+1) At] = Fj (At) F^(nAt) , F..(0) = I , j = 1,2,3,. 

G [ (n+1) At] = F 1 (At)G 1 (nAt) + G 1 (At)F 2 (nAt) , G^O) = 0 


G [(n+1)At] = F (At)G (nAt) + G (At)F (nAt) , G (0) = 0 


from which the integrals of Eq. (5-25) and (5-26) follow as 


/ e'^e^dt - F 2 [(n+1)At] Gjtn+DAt] (5-27) 


(n+1)At rT 

/ e A T D 3 e nT d T » F^[(n+1)At] G J(n+1)At.] (5-28) 


where At = (tf - t Q )/m and m is the total number of discrete time 


5.2.4.2 Vector Recursion Relationships for r[(n+1)At] 

Substituting Eq. (5-27) and (5-28) into Eq. (5-21), yields 


(n t 1)At -At T 
/ e 0 (t)C QCx*(t)dT 


F 3 [(n+1)At] G 1 [(n+1)At]s Q 


D 2 F^[(a+))Ac]G 2 [(n+1)At]s 
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= F^tn+DAt^F^AtjG^nAt) + G 1 (At)F 2 (nAt> }s q 

+ D 3 F^[(n+1)At]{F 3 (At)G 2 (nAt) + G 2 (At)F^nAt) }s q 

where s T = [(x*) T (X*) T ] 

o L o o J 

The number of mathematical operations required in the equation 
above can be minimized by defining the following three vector functions 
and associated recursion relationships: 

jc' [(n+1)At] = G 1 [(n+1) At]s Q = F 1 (At)x 1 (nAt) + G 1 ( At)x. (nAt) ; x 1 (0) = 0 

x 2 [(n+1)At] = G 2 [(n+1)At]s Q = F 3 <At)x 2 (nAt) + G 2 (At)x 3 (nAt); x 2 (0) = 0 

x 3 [(n+1)At] ■ F 2 [(n+1)At] s q = F 2 <At)x 3 (nAt); x 3 (0) = s q 

As a result, the recursive expression for the integral in Eq. 
(5-19), can be written as 

(n+1) At - 

/ e“ AT 0(T)C T QCxMr)dT = F^ta+OAt] xJtn+DAt] 
o 

+ D 3 F^[(n+1)At] x 2 ( (n+1) At ] 

where the recursion relationships for the matrix partitions G-|[(n+1)At] 
and G 2 [(n+ 1 )At] have been replaced by the vector recursion 
relationships for x^[(n+1)At], X 2 [(n+ 1 )At], and X 3 [(n+ 1 )At], 
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5.2.4.3 Recursion Relationship for 0[(n+1)At] 


* Subject to the recursion relationships above, the time histories 
for P and £ in Eq. (5-12) easily follow, where o[(n+1)&t] in Eq, (5-15) 
is given by 


0[ (n+1) At] = 0 ss + F 1 (At) [©(nAt) - © ss ]^<At) 

and 

e<o> - s M t »J(t f )[e £ - e 8 .]r 3 (t £ ) 

5.2.5 Illustrative Examples 

The specific model considered in this section (see Figs. 5-1 and 
5-2) consists of a rigid hub with four identical elastic appendages 
attached symmetrically about the central hub. In particular, the 
following idealizations are considered: (i) single-axis maneuvers? (ii) 
in-plane motion? (iii) anti-symmetric deformations? (iv) small linear 
flexural deformations? (v) only the linear time-invariant form of the 
equation of motion is considered? and (vi) the control actuators are 
modelled as concentrated torque generating devices. The distributed 
control system for the vehicle consists of: (i) a single controller in 
the rigid part of the structure, and (ii) a single controller in each 
elastic appendage. In the optimal control performance index, however, 
the control weighting matrices are adjusted in order to have the rigid 
body control provide the primary torque for maneuvering the vehicle, 
while the appendage controllers act principally as vibrational 
suppressors. 

For all cases, the following configuration parameters are assumed: 
the moment of inertia of f>>e undeformed structure about the spin axis, I, 
is 6764 kg-m 2 ? the mass/length ratio of the four identical elastic 
appendages, p, is 0.04096 kg/m? the length of each cantilevered 
appendage, L, is 35 m? the flexural rigidity of each cantilevered 
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appendage, El, is 1500 kg-m^/s 2 ; and the radius of the rigid hub, r, 
is 1 m. For simplicity, each appendage is assumed to have one 'ontroller 
located half-way along its span. In the integrations over the mass and 
stiffness distributions, the radius of the hub is not neglected in 
comparison to the appendage length. In the structural modeling 
equations, the following comparison functions have been adopted as 
"assumed modes" 

<(> p (x-r) « 1 - cos[pir(x-r)/L] + 0.5(-1 ) P+1 [pir(x-r)/L] 2 

(p = 1,2,...,n) (5-29) 

which satisfy the geometric and physical boundary conditions. 

4> i = 6' I = A* ' I T = d>* ' ' I T = 0 (5-30) 

T p|x=r T p|x=r T p |x=>r+L T p |x=r+L 

of a clamped free appendage, where (*)' = d(*)/dx. In addition, full 
state feedback is assumed for the results of this section. 

The output vector is assumed to be given by 



where 

y^ a central hub angular velocity 



y^ = four times the angular velocity at each apperdage controller 
location relative to the central hub 

y^ a hub angular position 

y 4 = four times the tip deflection of each appendage 
y*. = controls and control-rates 
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The number of elements in y 2 is given by the number of controllers on 
each of the four appendages. The output sub-vector y,. exists when 
control-rates are penalized. The elements of y g correspond to the 
additionally augmented states. 

Referring to’Tables 5-1 a and 5-2, the graphical summaries of the 
states and controls are discussed qualitatively in what follows. 

Case 1 (Fig. 5-3) presents a terminal controller example maneuver 
in which two flexible modes, are controlled in addition to the rigid body 
rotation. There are six states in this problem, all of which are 
specified at the final time. The time histories for this case are 
virtually identical to the results obtained from an open-loop controller 
using the same weight matrices in the performance index. The jump 
discontinuities in the control torques at the initial and final times are 
characteristic of controllers with no penalties on the control-rates in 
the performance index. As shown in Reference 5-1, the presence of these 
discontinuities tends to excite both modeled and unmodeled high frequency 
structural modes. 

Case 2 (Fig. 5-4) presents the same maneuver as in Case 1, except 
that the performance index now includes penalties on the first and second 
time derivatives of the control torques [Ref. 5-2). Including these 
penalties allows the control torques and torque-rates to be specified at 
the initial and final times, thereby eliminating the terminal jump 
discontinuities in the control profiles, uc this particular control 
problem, the states are augmented to include the control and first time- 
derivative of the control, while the second time-derivative of the 
control becomes the commanded input to the augmented system. On 
comparing with the results of Case 1, one finds much smoother modal 
amplitude and control torque time histories, although the peak torques 
and tip deflection are somewhat higher than those in Case 1. 
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Case 3 (Fig. 5-5) presents the results of a spin-to-rest maneuver 
in which the weight matrices and final conditions are identical with 
those of Case 2. The initial conditions are different from those in Case 
2, but the time-varying feedback gains are identical. As the plots in 
Figure 5-7 show, the terminal controller has little difficulty in 
bringing the system to the required final conditions. 

Case 4 (Fig. 5-6) presents the results of a spin-to-rest maneuver 
where off-nominal plant parameters are used. In particular, the 
moment of inertias for the hub and appendages are increased by 10% and 
the appendage mass per unit length is decreased by 10%, while the 
time-varying feedback gains are computed using the nominal structural 
parameters. Although the structural parameters are changed by only 10%, 
the overall effect, on the entire structure is quite large. Specifically, 
the first and second mode eigenvalues, u> 2 , are increased by 16% and 22% 
respectively, while the ratio of the hub moment of inertia to appendage 
moment of inertia is increased by 22%. Nevertheless, the plots show that 
most of the state variables reach their prescribed terminal boundary 
conditions. 

In Cases 1 through 4 the results have been compared with the 
corresponding open-loop maneuvers, and have been found to be identical 
within plotting accuracy, as expected. Furthermore, the closed form 
solutions for P(t) and $(t) have also been compared with the backward 
integrated solutions of Eqs. (5-10) and (5-11) for verification. 

5.3 Closed Form Solution for the Terminal Tracking Problem 

In this subsection the problem of maneuvering a flexible 
spacecraft through a large angle is considered, whsre the terminal state 
is subject to constraints. The basic approach differs from the standard 
linear quadratic regulator design, in that a quadratic penalty term on 


5-14 






the final value of the state is not included in the performance index. 

As a result, the numerical difficulties associated with handling large 
terminal penalty matrices are not encountered. 

The optimal control problem of this section is specified by 
defining a performance index which consists of an integral of quadratic 
forms in the state, control, and control-rates, where the terminal state 
is subject to constraints. In particular, the' terminal constraints 
replace the conventional weighted quadratic form in the terminal states. 
The necessary conditions for this problem lead to three nonlinear 
Riccati-like differential equations [Ref. 5-6]. One equation is the 
standard Riccati equation, which possesses a well-known solution in terms 
of a steady-state plus transient term (see Eq. (5-13)). On the other 
hand, the remaining two equations are coupled to the standard Riccati 
equation. In Sections 5.3.3 and 5.3.4, new closed form solutions are 
presented for the two auxiliary Riccati-like equations. 


5.3.1 Optimal Terminal Controller 

The optimal terminal control problem is formulated by finding the 
control inputs u(t) to minimize 

fc f 

J = J f (xVqfx + u T Ru) dt (5-32) 

t o 

for the system 

x = Ax + Bu , given x(t Q ) (5-33) 


with outputs 
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(5-34) 


y = Fx 

where x, u. A, B, Q, and R are defined following Eq. (5-3), F is the 
measurement influence matrix, and the state is subject to the specified 
terminal constraints 


x <t ) = X i , i = !,•••,q (q £ n) (5-35) 


Of particular interest is the fact that the performance index of. Eq. 
(5-32) does not contain a terminal weight matrix which penalizes the 
final values of the state. 


5.3.1.1 Necessary Conditions for the Terminal Tracking Problem 

As shown in Reference 5-6, the necessary conditions defining the 
optimal solution are given by the following coupled Riccati-like matrix 
differential equations 


where iji = (x 1 , 
given by 


A T P - PBR -1 B T P + F T QF = 0 ; P(t f ) = 0 (5-36) 

-PBR“ 1 B T )S = 0 ; S T (t f ) » (3<|)/3x)| t (5-37) 

G » S T BR -1 B T S ; G(t f ) = 0 (5-38) 

x_, •••, x ] and the optimal continuous feedback law is 
2 q 


u(t) = -C(t)x(t) - D(t)<|) 


(5-39) 


where 
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(5-40) 


(5-41) 


The control vector u of Eq. (5-39) takes the state vector from x(t Q ) at 
time t Q to at time tf, while minimizing the cost function of Eq. 

(5-32) . 

In the next section closed form solutions are presented for Eqs. 
(5-36), (5-37), and (5-38), thus reducing the solution for the feedback 
gains to the direct computation of algebraic equations without numerical 
integration. 

5,3,2 Closed Form Solution for the Time-Varying Riccati Equation 

P(t) » P^+e^U) (5-42) 

ss 


where P gs is the solution for the algebraic Riccati equation 


ss ss 


and the solution for 0(t) follows as 


A(t-t ) 

9 (t) = 9„ -. ie s3 * p;;>. 




(.) _ -IT 

where e is the exponential matrix, A = A - BR B P gg , and the solution 

for 0 gs is defined by Eq. (5-16). 


5.3.3 Closed Form Solution for S(t) 

The new solution for S(t) in Eq. (5-37) follows on assuming the 
product form solution 
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S(t) 


0 -1 (t)S (t) 


(5-45) 


where 0 is defined by Eq. (5-44) and S c is to be determined. 

Substitution of Eq. (5-45)- into Eq. (5-37) leads to the following 
linear constant coefficient matrix differential equation for S c (t): 

S (t) - A S (t) * 0; S (t.) = -P" 1 S(tJ (5-46) 

c c c f ss f 

from which it follows that the solution for S(t) in Eq. (5-45) can be 
written as: 


A(t-t^) 

S(t) =• -0 ' (t)e 1 P S(t_) 

ss f 

• - -T -IT 

where the differential equation for 0 is given by: 0 = A0 + 0A •- BR B . 
The solution for S c (t) in Eq. (5-46) can be shown to be 

A(t-t ) 

S (t) = -e p"' S(tJ (5-47) 

c ss f 

5.3.4 Closed Form Solution for G(t) 

The solution for G(t) in Eq. (5-38) can be shown to be: 

G(t) = S T (t)0(t)S(t) + S T (t,)P -1 S(t_) (5-48) 

f ss f 

which can be easily verified by direct differentiation. 

5.3.5 Calculation of the Optimal Control 

The control is computed by writing Eq. (5-39) in the form 


5-18 






u(t) 


= -C(t) x(t) - D(t)i|» 

- -r" 1 B T [p x(t) + e" 1 (t)x(t) + S(t)xft)] (5-49) 

1 ss 1 1 

where x-|(t) is obtained by solving the following linear equation: 

GCtJx^t) = - S(t)x-(t) 


5.3t6 Example Maneuvers 

The vehicle configuration and system parameters of Section 5.2.5 
are assumed for the example maneuvers of this section. Furthermore, full 
state feedback is assumed for the results of this section and the output 
vector is assumed to be defined by Eq. (5-31). 

Referring to Tables 5-1b and 5-3, the graphical summaries of the 
states and controls are discussed qualitatively in what follows. 

Case 1 (Fig. 5-7) presents a terminal controller example maneuver 
in which two flexible modes are controlled in addition to the rigid bod.' 
rotation. There are six states in this problem, all of which are 
specified at the final time. As a result, P, S, and G in Eqs. (5-42), 
(5-45), and (5-48) are 6x6 matrices. The time histories for this case 
are virtually identical to the time histories obtained from the 
corresponding open-loop controller, using the same weight matrices in the 
performance index. The observed jump discontinuities in the control 
torques at the initial and final times are characteristic of this type of 
control design. 

Case 2 (Fig. 5-8) presents a spin-up maneuver with two flexible 
modes controlled, in which the final angle is free to be determined by 
the controller. As a result, the P, S, and G matrices have the following 
dimensions: P(6x6), S(6x5), and G(5x5). In particular, the final angle 
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selected is exactly the angle selected by the analogous open-loop free 
final angle transversality conditions [Ref. 5-1]. 

Case 3 (Fig. 5-9) presents the same maneuver as in Case 1, except 
that the performance index now includes penalties on the first and second 
time derivatives of the control torques (Ref. 5-2]. Including these 
penalties allows the control torques and torque-rates to be specified at 
the initial and final times, thereby eliminating the terminal jump 
discontinuities in the control profiles. For this particular control 
problem, the states are augmented to include the control and first time- 
derivative of the control, while the second time-derivative of the 
control becomes the commanded input to the augmented system. On 
comparing with the results of Case 1, one finds a much smoother modal 
amplitude and control torque time history, although the peak torques and 
tip deflection are somewhat higher than those in Case 1. Hie dimensions 
of the P, S, and G matrices are all 10x10 for Case 3. 

Case 4 (Fig. 5-10) presents the same maneuver as in Case 2, with 
the additional penalty on the control-rates. As in Case 2, all the 
states, controls, and control-rates are specified at the final time for 
the terminal controller, except for the maneuver angle, which is not 
specified at the final time. As a result, the P, S, and G matrices have 
the following dimensions: P(10x10), S(10x9), andG(9x9). On comparing 
the results with Case 2, the modal amplitude and control torque time 
hiscories are found to be much smoother, although the peak values are 
higher. Moreover, the terminal jump discontinuities in the control 
torque shown in Case 2 are eliminated in Case 4, thus making it less 
susceptible to control spillover. 

Case 5 (Fig. 5-11) presents the results of a maneuver in which the 
weight matrices and final conditions are identical with those of Case 3. 
The initial conditions are different from those in Case 3, but the 
time-varying feedback gains are identical. As the plots in Figure 5-10 
show, the terminal controller has little difficulty in bringing the 
system to the required final conditions. 
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Case 6 (Fig. 5-12) presents the results of a maneuver in which the 
moment of inertias of the hub and appendages ''re increased by 10% and the 
appendage mass per unit length is decreased by 10%, while the 
time-varying feedba'ck gains are computed using the nominal structural 
parameters. Although the structural parameters are changed by only 10%, 
the overall effect or. the entire structure is quite large. Specifically, 
the first and second mode eigenvalues, to 2 , are increased by 16% and 
22%, respectively, while the ratio of the hub moment of inertia to 
appendage moment of inertia is increased by 22%. However, the plots show 
that with the terminal controller, most of the state variables reach 
their prescribed terminal boundary conditions. 

The results of in Cases 1 through 5 have been compared with the 
corresponding open-loop maneuvers, and have been found to be identical 
within plotting accuracy, as expected. Furthermore, the closed form 
solutions for P(t), S(t), and G(t) have also been compared with the 
backward integrated solutions for verification. 


5.4 An Analytic Solution for the State Trajectories of a Feedback 

Control System 

As part of the normal process of a control system design, the 
analyst typically is interested in determining the state trajectories for 
the controlled system. In practice, this process is straightforward, 
since the feedback form of the control can be introduced in the equation 
of motion and numerically integrated. Nevertheless, this process can be 
computationally intensive if either time-varying control gains are used 
or small integration step sizes are required by the presence of high 
frequency system dynamics. 

In an effort to overcome the computational difficulties listed 
above, a change of variables is presented in this subsection for the 
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standard closed-loop system dynamics equation, which permits a closed 
form expression to be obtained for r'ne state trajectories. 


5.4.1 Optimal Control Problem 

The fixed time linear optimal control problem is formulated by 
finding the control inputs u(t) to minimize 

t f 

J = ~ x^F T SFx f + J (x T F T QFx + u T Ru) dt (5-50) 

fc o 

for the system 

x = Ax + Bu , given x(t Q ) (5-51) 

y = Fx (5-52) 

where x, u, A, B, F, Q, R, and S are defined following dg. (5-3). 

As shown in Reference 5-6, the optimal control Is. given by 

u (t) a -R"°B T P(t)x(t) (5-53) 

where P is the solution for the differential matrix Riccaci equation 

P = -AP - PA T + PBR -1 B T P - Q j P(t f ) = S (5-54) 

Upon introducing Eq. (5-53) into Eq. (5-51), the standard closed- 
loop system dynamics equation follows as 
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x(t) = [a - BR -1 B T P(t)lx(t) ; x = x(t ) (S™55) 

1 J o o 

(It should be observed that the equation above is the adjoint to the 
homogeneous part of Eq. (5-11).) 

To obtain the solution for Eq. (5-55), the following closed form 
solution for P(t) [Refs. 5-1, 5-3] is introduced in Eq. (5-55): 

P(t) = P ss +0 " 1(t) (5 " 56) 

where P ss is is defined by Eq. (5-43). The variable 0(t) is given by 

*<»V, t -i , » T <'-y 

S<«> - S ss *. (CSf-P^ -Sje (5-57) 

- -IT (•) 

where A = A - BR BP ,e is the exponential matrix, and the solution 
ss 

for 0 gs is defined by Eq. (5-16). 

Substituting Eq. (5-56) into Eq. (5-55), yields the modified form 
of the closed-loop system dynamics equation 

x(t) = [A - BR“ 1 B T 0” 1 (t)]x(t) / x = x(t ) (5-58) 

v i o O 

where it follows that the equation above is nonautonomous. 


5.4.2 Change of Variables for x(t) 


To simplify Eq. (5-58) the following coordinate transformation for 
the dependent variable x(t) is introduced: 








x(t) 


0(t)r(t) 


(5-59) 


where 0(t) is given by Eq. (5-57) and r(t) is a vector function which is 
to be determined. 

Upon differentiating Eq. (5-59) we find 


(5-CO) 


or 

X = (A0 + 0A T - BR _1 B T )r + 0r (5-61) 

• ' • - -T -1 

where 0 in Eq. (5-60) has been replaced by: 0 = A0 + 0A - BR B. 

The differential equation for r is obtained by introducing Eqs. 
(5-59) and (5-61) into Eq. (5-58), leading to 

0(r + A T r) = 0 (5-62) 

from which it follows that the linear constant coefficient vector 
differential equation for r is given by 


(5-63) 


The solution for r follows as 


r(t) 


(5-64) 


Substituting Eq. (5-64) into Eq. (5-59) produces the desired 
solution for the state trajectories as 
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-T 

“Alt—t ) 


x(t) = 0(t)e r Q 


(5-65) 


5.4.3 Recursion Relationship for Evaluating the State at Discrete Times 

If the solution for x(t) is required at the discrete times 
=> t c + kAt (k » 1, • • »,N) for At = (tf - t c )/N, then Eq. (5-65) 
can be written as 

xft^) - ^r r + b k , k = O f ...,N (5-66) 


f _ TJ l" 1 _ 


A, » A. e 


-A* At 

k-1 e ' b k 


e“ AAt b k _ 1 ; k=1,««.,N 


5.4.4 Conclusions 

Closed form solutions for the feedback gains required by an 
optimal linear tracking controller and an optimal terminal controller 
have been developed. Results of example maneuvers have been shown which 
demonstrate the efficiency and validity of the formulations described. 

The use of control-rate penalties in each case have been shown to improve 
the overall system response. 
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In addition, a straightforward algorithm has been presented for 
generating the state trajectories for a feedback control system. The 
algorithm is computationally efficient in that no numerical integration 
is required and simple recursion relationships generate the desired 
solution at discrete times. Furthermore, this algorithm has significant 
potential if used in conjunction with algorithms which attempt to enhance 
system robustness, by iteratively refining the weighting matrices 
appearing in Eq. (5-50). 
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Table 

5-la. 

Test case 

linear tr« 

maneuver descriptions: 

aeking problem^ . 

No. of 

V fc o 


0 O 

6 f 

9 f 

Modes 

(sec) 

(rad) 

(rad/s) 

(rad) 

(rad/s) 

2 

15 

0 0 

0 

0.4 

0 

2 

15 

2 0 

0 

0.4 

0 

2 

15 

2 0.15 

-0.02 

0.4 

0 

2 

1.5 

2 0 

0 

0.4 

0 

Table 

5-1 b. 

Test case 

maneuver 

descriptions: 



terminal ' 

bracking ] 

ffoblem 



Case 

No. 

No. of 
Modes 

• V fc o 

(sec) 

k' 

0 o 

(rad) 

0 o 6 f 

(rad/s) (rad) 

0 f 

(rad/s) 

CPU 

(sec) 

1 

2 

15 

0 

0 

0 0.4 

0 

34 

2 

2 

15 

0 

0 

0 (free) 

0.04 

32 

3 

2 

15 

2 

0 

0 0.4 

0 

83 

4 

2 

15 

2 

0 

0 (free) 

0.04 

77 

5 

2 

15 

2 

0.15 

-0.02 0.4 

0 

83 

6 

2 

15 

2 

0 

0 0.4 

0 

83 

Notes: 








U 

k denotes 
which is 

the order of 
penalized in 

the highest time derivative of 
the performance index. 

the cont 

2. 

The CPU time is obtained using 
steps for each integration. 

an Amdahl 470-V/8 with 

150 time 































Weighting matrices for example maneuvers: 
terminal tracking problem. 
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Figure 5-5. Case 3, spin-to-rest maneuver, 2 modes, 5 controls, 
control-rate penalty with k * 2. 
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Case 4, rest-to-rest maneuver, off-nominal structural 
parameters, 2 modes, 5 controls, k * 2. 
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